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Abstract 

One of the most surprising consequences of quantum mechanics is the nonlocal correlation 
of a multi-particle system observable in joint-detection of distant particle-detectors. Ghost 
imaging is one of such phenomena. Taking a photograph of an object, traditionally, we need 
to face a camera to the object. But with ghost imaging, we can image the object by pointing 
a CCD camera towards the light source, rather than towards the object. Ghost imaging is 
reproduced at quantum level by a non-factorizable point-to-point image-forming correlation 
between two photons. Two types of ghost imaging have been experimentally demonstrated 
since 1995. Type-one ghost imaging uses entangled photon pairs as the light source. The non- 
factorizable image-forming correlation is the result of a nonlocal constructive-destructive inter- 
ference among a large number of biphoton amplitudes, a nonclassical entity corresponding to 
different yet indistinguishable alternative ways for the photon pair to produce a joint-detction 
event between distant photodetectors. Type-two ghost imaging uses chaotic light. The type- 
two non-factorizable image-forming correlation is caused by the superposition between paired 
two-photon amplitudes, or the symmetrized effective two-photon wavefunction, corresponding 
to two different yet indistinguishable alternative ways of triggering a join-detection event by 
two independent photons. The multi-photon interference nature of ghost imaging determines 
its peculiar features: (1) it is nonlocal; (2) its imaging resolution differs from that of classical; 
and (3) the type-two ghost image is turbulence-freeQ Ghost imaging has attracted a great deal 
of attention, perhaps due to these features for certain applications. Achieving these features, 
the realization of nonlocal multi-photon interference is a necessary condition. Classical sim- 
ulations, such as the man-made factorizable speckle-speckle correlation, can never have such 
features. 



1 Introduction 

Assuming an object that is cither self-luminous or externally illuminated, imagining each point 
on the object surface as a point radiation sub-source, each point sub-source will emit spherical waves 
to all possible directions. How much chance do we expect to have a spherical wave collapsing into 
a point or a "speckle" by free propagation? Obviously, the chance is zero unless an imaging system 
is applied. The concept of optical imaging was well developed in classical optics for this purpose. 
Figure [T] schematically illustrates a standard imaging setup. In this setup an object is illuminated 
by a radiation source, an imaging lens is used to focus the scattered and reflected light from the 
object onto an image plane which is defined by the "Gaussian thin lens equation" 
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where s a is the distance between the object and the imaging lens, Sj the distance between the 
imaging lens and the image plane, and / the focal length of the imaging lens. Basically this 
equation defines a point-to-point relationship between the object plane and the image plane: any 
radiation starting from a point on the object plane will "collapse" to a unique point on the image 
plane. It is not difficult to see from Fig. [I] that the point-to-point relationship is the result of 
constructive-destructive interference. The radiation fields coming from a point on the object plane 
will experience equal distance propagation to superpose constructively at one unique point on 
the image plane, and experience unequal distance propagations to superpose destructively at all 



1 For instance, any fluctuation of the refraction index or phase disturbance in the optical path has no influence 
to the type-two ghost image. 
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Figure 1: Optical imaging: a lens produces an image of an object in the plane denned by the 
Gaussian thin-lens equation 1/sj + l/s a = 1//. Image formation is based on a point-to-point 
relationship between the object plane and the image plane. All radiations emitted from a 
point on the object plane will "collapse" to a unique point on the image plane. 



other points on the image plane. The use of the imaging lens makes this constructive-destructive 
interference possible. 

A perfect point-to-point image-forming relationship between the object and image planes pro- 
duces a perfect image. The observed image is a reproduction, either magnified or demagnificd, 
of the illuminated object, mathematically corresponding to a convolution between the object dis- 
tribution function |A(p )| 2 (aperture function) and a ^-function which characterizes the perfect 
point-to-point relationship between the object and image planes: 

I(pi)= [ dp \A(p )\ 2 S(p +^) (2) 

J obj 171 

where I(f>i) is the intensity in the image plane, p Q and pi are 2-D vectors of the transverse coor- 
dinates in the object and image planes, respectively, and m — Si/s is the image magnification 
factor. 

In reality, limited by the finite size of the imaging system, we may never obtain a perfect point- 
to-point correspondence. The incomplete constructive-destructive interference turns the point-to- 
point correspondence into a point-to- "spot" relationship. The <5-function in the convolution of 
Eq. ([2]) will be replaced by a point-spread function: 

I{pi) = / dp \A(p )\ somb 2 [ \p a + — |], (3) 

Jobj s Q c m 

where the sombrero-like function, or the Airy disk, is defined as 

u \ 2Jl ^ 

somb(x) = , 

x 

and Ji(x) is the first-order Bessel function, and R the radius of the imaging lens, and R/s is 
known as the numerical aperture of the imaging system. The sombrero-like point-spread function, 
or the Airy disk, defines the spot size on the image plane that is produced by the radiation coming 
from point p a . It is clear from Eq. ^ that a larger imaging lens and shorter wavelength will result 
in a narrower point-spread function, and thus a higher spatial resolution of the image. The finite 
size of the spot determines the spatial resolution of the imaging system. 

Type-one and type- two ghost imaging, in certain aspects, exhibit a similar point-to-point 
imaging-forming function as that of classical except the ghost image is reproducible only in the 
joint-detection between two independent photodetectors, and the point-to-point imaging-forming 
function is in the form of second-order correlation, 

RM) = f dp \A(p )\ 2 G^(p ,p t ), (4) 

J obj 
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where Ri 2 (pi) is the joint-detection counting rate between photodetectors D\ and D 2 . Mathe- 
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matically, the convolution is taken between the aperture function of the object L4(p ) and a 
nontrivial poin-to-point second-order correlation function G^ 2 \p a , pi), corresponding to the prob- 
ability of observing a joint photo-detection event at coordinates p and pi. It is the special physics 
behind G^ 2 \p ,pi) made ghost imaging so special. 
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Figure 2: Schematic set-up of the first "ghost" image experiment. The experimental demon- 
strations of ghost imaging and ghost interference jl] in 1995 together stimulated the foundation 
of quantum imaging in terms of geometrical and physical optics. 



The first type-one ghost imaging experiment was demonstrated by Pittman et at in 1995 
[1] enlightened by the theoretical work of Klyshko [2]. The schematic setup of the experiment is 
shown in Fig. [2j A continuous wave (CW) laser is used to pump a nonlinear crystal to produce an 
entangled pair of orthogonally polarized signal (e-ray of the crystal) and idler (o-ray of the crystal) 
photons in the nonlinear optical process of spontaneous parametric down-conversion (SPDC). The 
pair emerges from the crystal collinearly with u> s = Ui = lo p /2 (degenerate SPDC). The pump is 
then separated from the signal-idler pair by a dispersion prism, and the signal and idler are sent 
in different directions by a polarization beam splitting Thompson prism. The signal photon passes 
through a convex lens of 400mm focal length and illuminates a chosen aperture (mask). As an 
example, one of the demonstrations used the letters "UMBC" for the object mask. Behind the 
aperture is the "bucket" detector package D\, which is made by an avalanche photodiode placed 
at the focus of a short focal length collection lens. During the experiment D\ is kept in a fixed 
position. The idler photon is captured by detector package D 2 , which consists of an optical fiber 
coupled to another avalanche photodiode. The input tip of the fiber is scannable in the transverse 
plane by two step motors (along orthogonal directions). The output pulses of D\ and D 2 , both 
operate in the photon counting mode, are independently counted as the counting rate of D± and 
D 2l respectively, and simultaneously, sent to a coincidence circuit for counting the joint-detection 
events of the pair. The single detector counting rates of D\ and D 2 are both monitored to be 
constants during the measurement. Surprisingly, a ghost image of the chosen aperture is observed 
in coincidences during the scanning of the fiber tip, when the following two experimental conditions 
are satisfactory: (1) D± and D 2 always measure a pair; (2) the distances s , which is the optical 
distance between the aperture to the lens, Sj, which is the optical distance from the imaging lens 
going backward along the signal photon path to the two-photon source of SPDC then going forward 
along the idler photon path to the fiber tip, and the focal length of the imaging lens / satisfy the 
Gaussian thin lens equation of Eq. ([IJ . 

Figure [3] shows a typical measured ghost image. It is interesting to note that while the size 
of the "UMBC" aperture inserted in the signal path is only about 3.5mmx7mm, the observed 
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Figure 3: Upper: A reproduction of the actual aperture "UMBC" placed in the signal beam. 
Lower: The image of "UMBC" : coincidence counts as a function of the fiber tip's transverse 
coordinates in the image plane. The step size is 0.25mm. The image shown is a "slice" at the 
half maximum value. 

image measures 7mm x 14mm. The image is therefore magnified by a factor of 2 which equals the 
expected magnification m = Si/s a . In this measurement s = 600mm and Sj = 1200mm. When 
Z?2 was scanned on transverse planes other than the ghost image plane the images blurred out. 

The experiment was immediately given the name "ghost imaging" by the physics community 
due to its nonlocal feature. In the language of Einstein-Podolsky- Rosen (EPR) [3], the non- 
factorizabk[^] point-to-point image-forming correlation 

G (2) (/W») ~ S(p +pi/m) (5) 

observed in this experiment represents a nonlocal behavior of a measured pair of photons: neither 
the signal photon nor idler photon "knows" precisely where to go when the pair is created at 
the source. However, if one of them is observed at a point on the object plane, the other one 
must arrive at a unique corresponding point on the image plane j^] Although questions regarding 
fundamental issues of quantum theory still exist, the experimental demonstration of ghost imaging 
[T] and ghost interference jl] in 1995 together stimulated the foundation of quantum imaging in 
terms of geometrical and physical optics. 

Type-two ghost imaging uses chaotic radiation sources. Different from type-one, the non- 
factorizable point-to-point image-forming correlation between the object and image planes is only 
partial with at least 50% constant background, 

G (2) (p ,Pi) ~\ + 8{p -pi). (6) 

The first near-field lensless ghost imaging experiment was demonstrated by Scarcelli et al in 2005 
and 2006 |S] [Hj after their experimental demonstration of two-photon interference of chaotic light 
in 2004 [7] . Figure [4] illustrates an improved setup of the type-two ghost imaging experiment by 
Meyers et al. [8]. The thermal radiation of a chaotic source, which has a fairly large size in the 
transverse dimension, is split into two by a 50% — 50% beamsplitter. One of the beams illuminates 

2 Statistically, a factorizable correlation function G^ 2 ' (ri , ti ; r-2, (2) = CrW (n, ti) C?W (r2, tz) characters inde- 
pendent radiations at space-time (n,ii) and (r2,*2)- In ghost imaging, the light on the object plane and the light 
at the CCD array is described by a non-factorizeable point-to-point image-forming function, indicating nontrivial 
statistical correlation between the two measured intensities. 

3 The ghost imaging experiment is thus considered a demonstration of the historical Einstein-Podolsky-Rosen 
(EPR) experiment. 
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Figure 4: Near-field lensless ghost imaging of chaotic light demonstrated by Meyers et al. D 2 
is a "bucket" photon counting detector that is used to collect and count all random scattered 
and reflected photons from the object. The joint-detection between D 2 and the CCD array 
is realized by a photon-counting-coincidence circuit. D 2 is fixed in space. The counting rate 
of D 2 and the un-gated output of the CCD are both monitored to be constants during the 
measurement. Surprisingly a 1:1 ghost image of the object is captured in joint-detection 
between D 2 and the CCD, when taking Z\ = z<i- The images "blurred out" when the CCD is 
moved away from z\ — z 2 , either in the direction of z\ > z 2 or z\ < z 2 . 




Figure 5: Ghost image of a toy soldier model. 



a toy soldier as shown in Fig. |1] The scattered and reflected photons from the solider (object) are 
collected and counted by a "bucket" detector D 2 . In the other beam a high resolution CCD array 
operating at the photon counting regime, is placed toward the radiation source for joint-detection 
with the "bucket" detector D 2 . The counting rate of D 2 and the un-gated output of the CCD are 
both monitored to be constants during the measurement. Surprisingly, a 1:1 ghost image of the 
toy soldier is captured in the joint-detection between D 2 and the CCD, when taking z\ — z 2 . The 
1:1 ghost image of the toy soldier is shown in Fig. [5| The images "blurred out" when the CCD is 
moved away from z\ = z 2 , either to the side of z\ > z 2 or z\ < z 2 . 

There is no doubt that chaotic radiations propagate to any transverse plane in a random 
and chaotic manner. A brief discussion for Fresnel free-propagation is given in the appendix. In 
the lensless ghost imaging experiment, a large transverse sized chaotic: light source, as shown in 
Fig. [4] is usually used for achieving better spatial resolution. The source consists a large number 
of independent point sub-sources randomly distributed on the source plane. Each point sub-source 
may randomly radiate independent spherical waves to the object and image planes. Due to the 
chaotic nature of the source there is no interference between these sub-fields. These independent 
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sub-intcnsitics simply add together, yielding a constant total intensity in space and in time on 
any transverse plane. In the lcnslcss ghost imaging setup, there is no lens applied to force these 
spherical waves collapsing to a point or a "speckle" , and there is no chance to have two identical 
copies of any "speckle" of the source onto the object and image planes. What is the physical 
cause of the point-to-point image-forming correlation? Although the non-factorizable point-to- 
point correlation between the object and image planes is only partial, the type-two ghost imaging 
looks more surprising than type-one because of the nature of the light source. Unlike the signal- 
idler photon pair, the jointly measured photons in type-two ghost imaging are just two independent 
photons that fall into the coincidence time window by chance only. Nevertheless, analogous to EPR, 
the non-factorizable partial point-to-point correlation represents a nonlocal behavior of a measured 
pair of independent photons: neither photon-one nor photon-two "knows" precisely where to go 
when they are created at each independent sub-sources; however, if one of them is observed at 
a point on the object plane, the other one has twice greater probability of arriving at a unique 
corresponding point on the image plane j^] 

We have concluded and will show that the partial point-to-point correlation between the object 
and image planes in type-two ghost imaging is the result of two-photon interference. Similar to that 
of type-one, it involves the nonlocal superposition of two-photon amplitudes, a nonclassical entity 
corresponding to different yet indistinguishable alternative ways of triggering a joint-detection event 
[§]. Different from that of type-one, the joint-detection events observed in type-two ghost imaging 
are triggered by two randomly distributed independent photons. It is interesting to see that the 
quantum mechanical concept of two-photon interference is applicable to "classical" thermal light 
In fact, this is not the first time in the history of physics we apply quantum mechanical concepts 
to thermal light. We should not forget Planck's theory of blackbody radiation originated the 
quantum physics. The radiation Planck dealt with was thermal radiation. Although the concept 
of "two-photon interference" comes from the study of entangled biphoton states [9], the concept 
should not be restricted to entangled systems. The concept is generally true and applicable to any 
radiation, including "classical" thermal light. The partial point-to-point correlation of thermal 
radiation is not a new discovery either. The first set of temporal and spatial far-field intensity- 
intensity correlations of thermal light was demonstrated by Hanbury Brown and Twiss (HBT) 
in 1956 |10||11]. The HBT experiment created quite a surprise in the physics community and 
lead to a debate about the classical or quantum nature of the phenomenon [H][12]. Although 
the discovery of HBT initiated a number of key concepts of modern quantum optics, the HBT 
phenomenon itself was finally interpreted as statistical correlation of intensity fluctuations and 
considered as a classical effect. It is then reasonable to ask: Is the near-field type-two ghost 
imaging with thermal light a simple classical effect similar to that of HBT? Is it possible that 
the ghost imaging phenomenon itself, including the type-one ghost imaging of 1995, is merely a 
simple classical effect of intensity fluctuation correlation? [13] [14] [15] [16] This article will address 
these important questions and explore the multi-photon interference nature of ghost imaging. 

To explore the two-photon interference nature, we will analyze the physics of type-one and 
type-two ghost imaging in five steps. (1) Review the physics of coherent and incoherent light prop- 
agation; (2) Review classical imaging as the result of constructive-destructive interference among 
electromagnetic waves; (3) analyze type-one ghost imaging in terms of constructive-destructive 
interference between the biphoton amplitudes of an entangled photon-pair; (4) analyze type-two 
ghost imaging in terms of two-photon interference between chaotic sub-fields; and (5) discuss the 
physics of the phenomenon: whether it is a quantum interference or a classical intensity fluctuation 
correlation. 

4 Similar to the HBT correlation, the contrast of the near-field partial point-to-point image-forming function is 
50%, i.e., two to one ratio between the maximum value and the constant background, see Eq. ( |33[ ) ■ 

5 There exist a number of definitions for classical light and for quantum light. One of the commonly accepted 
definitions considers thermal light classical because its positive P-function. 
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2 Classical Imaging 

To understand the multi-photon interference nature of ghost imaging, it might be helpful to 
see the constructive-destructive interference nature of classical imaging first. We start from a 
typical classical imaging setup of Fig. [6] and ask a simple question: how does the radiation field 
propagate from the object plane to the image plane? In classical optics such propagation is usually 
described by an optical transfer function h(r — r , t — t ). We prefer to work with the single-mode 
propagator, namely the Green's function, g(k, r — 1*0, t — to) [H][IB], which propagates each mode 
of the radiation from space-time point (ro,£o) to space-time point (r,t). We treat the field E(r,i) 
as a superposition of these modes. A detailed discussion about <?(k,r — ro,t — to) is given in the 
Appendix. It is convenient to write the field E(r, t) as a superposition of its longitudinal and 
transverse modes under the Fresnel paraxial approximation, 

E(p,z,t) = J dnduj E(k,uj) g(K,uj; p, z) e~ luJt , (7) 

where E{H, u>) is the complex amplitude for the mode of frequency u> and transverse wave- vector k. 
In Eq. we have taken zq = and to = at the object plane as usual. To simplify the notation, 
we have assumed one polarization. 
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Figure 6: Typical imaging setup. A lens of finite size is used to produce a magnified or 
demagnified image of an object with limited spatial resolution. 

Based on the experimental setup of Fig. [6] and following the Appendix, g(n, ui; p, z) is found to 



g(K,u;pi,s + Si) 
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(8) 



where p 01 pi, and pi are two-dimensional vectors defined, respectively, on the object, lens, and 
image planes. The first curly bracket includes the aperture function A(p a ) of the object and 
the phase factor e tK ' p ° contributed at the object plane by each transverse mode k. The terms 
in the second and fourth curly brackets describe free-space Fresnel propagation-diffraction from 
the source/object plane to the imaging lens, and from the imaging lens to the detection plane, 
respectively. The Fresnel propagator includes a spherical wave function e l ^^ Zi ~ Zk ^ > /{zj — z^) and a 
Fresnel phase factor e tul ^ j ^^ 2 ^ 2c< - Zj ^ Zk K The third curly bracket adds the phase factor introduced 
by the imaging lens. 

We now rewrite Eq. d8l into the following form 



g(K,cj;pi,z = s Q + = 



(27TC) 2 S S, 



dp A(p ) e 



:\Po\ 2 e tK-Pa 



obj 



(9) 



lei is 
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The image plane is defined by the Gaussian thin- lens equation of Eq. Q. Hence, the second integral 
in Eq. ^ reduces to, for a finite sized lens of radius R, the so-called point-spread function, or the 
Airy disk, of the imaging system: 

dpt e-^fe+t^' = 2J ^1 = SO mb(x), (10) 

where the sombrero-like function somb(x) = 2Ji(x)/x with argument x — [— — \p D + pi/m\) has 
been defined in Eq. (j3j. Eq. (10 1 indicates a constructive interference. 



Substituting Eqs. (JsJ) and ( 10 1 into Eq. ^ enables one to obtain the classical self-correlation 
function of the field, or, equivalently, the intensity on the image plane 

I(pi,Zi,ti) = ( E*(p i ,z l , ti) E(pi, Zi,ti) ), (11) 

where (...) denotes an ensemble average. To simplify the mathematics, monochromatic light is 
assumed as usual. 

Case (I): Incoherent imaging. The ensemble average yields zeros except when k — k' . The 
image is thus 

I(pi) oc [ dp \A(p )\ 2 somb 2 [--\p +^\}. (12) 
J 1 s c m 

An incoherent image, magnified by a factor of m, is thus given by the convolution between the mod- 
ulus square of the object aperture function and the point-spread function. The spatial resolution 
of the image is determined by the finite width of the |som&| 2 -function. 

Case (II): Coherent imaging. The coherent superposition of the k modes in both E*(p>i, r) and 
E(pi,r) results in a wavepacket. The image, or the intensity distribution on the image plane, is 

I(pi)cx [ dp A(po)e^^ 2 S omb[--\po+^\} 2 . (13) 
Jobj s c m 

A coherent image, magnified by a factor of m, is thus given by the modulus square of the convolution 
between the object aperture function (multiplied by a Fresnel phase factor) and the point-spread 
function. 



For Si < s and s Q > /, both Eqs. (12 1 and (13 1 describe a real demagnified inverted image 



In both cases, a narrower somfe-function yields a higher spatial resolution. Therefore the use of 
a larger imaging lens and shorter wavelengths will improve the spatial resolution of an imaging 
system. 

3 Biphoton and type-one ghost imaging 

In this section we analyze type-one ghost imaging. Type-one ghost imaging uses entangled 
photon pairs such as the signal-idler biphoton pairs of SPDC [13] [5] • The nearly collinear signal- 
idler system generated by SPDC can be described, in the ideal case, by the following entangled 
biphoton state [§]: 

|*)=* J dK s dKiS(K s + Hi) J (L) s doJiS(uj s u)p)a^(K s ,u} s )a^(Ki,u!i)\0), (14) 
where toj, iij (j = s,i,p), are the frequency and transverse wavevector of the signal, idler, and 



pump, respectively. For simplicity a CW single mode pump with k p = is assumed. Eq. (14) 
indicates that the biphoton state of the signal-idler pair is an entangled state. The single-photon 
state of the signal and the idler can be evaluated by taking a partial trace of its twin, 



p s = tri | * )(* | = / dn s duj s at (k s , uj s )\ 0)(0 \o,(k s ,uj s ), 

p t = tr s |*)(*| = / dKiduJi a t (K l ,w l )|0)(0|a(K J ,cj l ). (15) 
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Although the signal-idler system is in a pure state, the state of the signal photon and the idler 
photon, respectively, are both mixed states. 

Let us imagine a measurement in which two point- like photon counting detectors (D\ and D2) 
are placed at the output plane of an SPDC source for the detection of the signal photon and the 
idler photon, respectively, and for the joint-detection of the signal-idler pair. The probability of 
observing a photo-detection event in the SPDC output plane pj at time tj, j = s, i, is calculated 
from the first-order photo-detection theory of Glauber [20 

& x) [pj,ti) = trpE^ipj^E^ip^tj), (16) 

where we have chosen zj = for the SPDC output plane as usual. It is easy to find that 

G^\p s ,t s ) — constant, G lyl \pi,ti) — constant, (17) 

which means that the signal photon and the idler photon both have equal probability to be observed 
at any position in the output plane of the SPDC at any time. The probability of observing a joint- 
detection event between D± and D2 located at p s and pi in the SPDC output plane of z s — Zi — 
is calculated from the second-order photo-detection theory of Glauber [2D] : 

G^{p sl t s -p l: U) 
= \(0\El + \p i ,t i )E^(p s ,t s )\^)\ 2 

= \^{p s ,t s ;pi,ti)\ 2 , (18) 

where ^(p s , t s ; Pi, U) is defined as the effective biphoton wavefunction. The transverse spatial part 
of the effective biphoton wavefunction is easily calculated to be: 

V{p s ,Pi)~5{p s -pi), (19) 



under the condition t s ~ tj. Equations (|14[), (17), and (19 1 suggest that the entangled signal- idler 



photon pair is characterized by the EPR correlation [3] in transverse momentum and transverse 
position; hence, similar to the original EPR state, we have [21 : 

A(k s + Zi) = & A(p s -pi) = i) (20) 
with Ak s ~ 00, A.Ki ~ 00, Ap s ~ 00, Api ~ 00. 

In EPR's language, the signal photon and the idler photon may come from any point in the output 
plane of the SPDC. However, if the signal (idler) is found in a certain position, the idler (signal) 
must be observed in the same position, with certainty (100%). Simultaneously, the signal photon 
and the idler photon may have any transverse momentum. However, if a certain value and direction 
of the transverse momentum of the signal (idler) is observed, the transverse momentum of the idler 
(signal) will be uniquely determined with equal value and opposite direction. 

The EPR (5-functions, S(p s — pi) and 5{k s + Kj) in transverse position and momentum, are the 
key to understanding the ghost imaging experiment of Pittman et al. of 1995. 6(p s — pi) indicates 
that the signal-idler pair is always emitted from the same point on the output plane of the biphoton 
source. Simultaneously, 6(k s + iti) defines the angular correlation of the pair: the signal-idler pair 
always exists at roughly equal but opposite angles relative to the pump for degenerate SPDC. 
This then allows for a simple explanation of the experiment in terms of "usual" geometrical optics 
in the following manner: we envision the nonlinear crystal as a "hinge point" and "unfold" the 
schematic of Fig. [2] into the Klyshko picture [5] of Fig. [7j The signal-idler biphoton amplitudes 
can then be represented by straight lines (but keep in mind the different propagation directions) 
and therefore the image is reproduced in coincidences when the aperture, lens, and fiber tip are 
located according to the Gaussian thin lens equation of Eq. ([T]). The image is exactly the same 
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Figure 7: An unfolded schematic of the 1995 ghost imaging experiment, which is helpful for 
understanding the physics. Since the biphoton "light" propagates along "straight lines" , it is 
obvious that any point on the object plane corresponds to a unique point on the image plane. 
Although the placement of the lens, the object, and detector D2 obeys the Gaussian thin lens 
equation, it is important to notice that the geometric rays in the figure actually represent the 
biphoton amplitudes of an entangled photon pair. The point-to-point correspondence is the 
result of a constructive-destructive interference of these biphoton amplitudes. 



as that one would observe on a screen placed at the fiber tip if detector Di were replaced by a 
point-like light source and the nonlinear crystal by a reflecting mirror. 

Comparing the "unfolded" schematic of the ghost imaging experiment with that of the classical 
imaging setup of Fig. [l] it is not difficult to find that any "light point" on the object plane has 
a unique corresponding "light point" on the image plane. This point-to-point correspondence is 
the result of the constructive-destructive interference among these biphoton amplitudes that are 
illustrated as the geometrical rays in Fig. [7j Similar to the situation in classical imaging, these 
biphoton amplitudes which experience equal optical path propagation will superpose constructively 
at each pair of one-to-one points of the object plane and the image plane for a joint-detection event, 
while these that experience unequal distance propagation will superpose destructively at all other 
points on the object and image planes. The use of the imaging lens makes this constructive- 
destructive interference possible. It is this unique point-to-point EPR correlation that makes the 
"ghost" image of the object-aperture function possible. Despite the completely different physics 
from classical optics, the remarkable feature is that the relationship between the focal length 
/ of the lens, the aperture's optical distance s , and the image's optical distance Sj, satisfies 
the Gaussian thin lens equation of Eq. ([I]). It is worth emphasizing again that the geometric 
rays in Fig. [7] represent the biphoton amplitudes of a signal-idler photon pair, and the point-to- 
point correspondence is the result of the constructive-destructive interference of these biphoton 
amplitudes. 

We now calculate G@\pi, fc) for the "ghost" imaging experiment in detail, where p\ and P2 
are the transverse coordinates of the point-like photodetector D\ and D2, on the object and image 
planes, respectively. We will show that there exists a (5-function-like point-to-point correlation 
between the object and image planes, 8{p\ — fy/m). We will then show how the object function 
of A(p ) is transferred to the image plane as a magnified image A(p2/m). 

We first calculate the effective biphoton wavefunction ^(pi, Zi,t\; p2, %2> £2)5 as defined in 
Eq. ( |18[ ). By inserting the field operators into ^{p\,Zi,ti;p2,Z2,t2), and considering the com- 
mutation relations of the field operators, the effective biphoton wavefunction is calculated to be 



$(pi,Zx,ti;p 2 ,Z2,t2) = * / dK s dKiS(fi s + Ki) j duj s dDiS(u a + uji - uj p ) (21) 

x g(% a ,u a ;pi,zi) e~^ stl g(Hi, uj 1 ; p 2 ,z 2 ) e" 1 ^ 2 . 
Equation (21 1 indicates a coherent superposition of all the biphoton amplitudes shown in Fig. [7] 
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Next, we follow the unfolded experimental setup of Fig. [8] to establish the Green's functions 
g(its) w s , pi, zi) and <?(/tj,Wj, p2, z%). In arm-1 the signal propagates freely over a distance d\ 
from the output plane of the source to the imaging lens, passes an object aperture at distance 
s , and then is focused onto photon-counting detector D\ by a collection lens. We will evaluate 
g(n s ,uj s , pi, Z\) by propagating the field from the output plane of the biphoton source to the object 
plane. In arm-2 the idler propagates freely over a distance di from the output plane of the biphoton 
source to a point-like detector £>2- gi^i, <^i, P21 z%) is thus a free propagator. 




Figure 8: In arm-1 the signal propagates freely over a distance d\ from the output plane of the 
source to the imaging lens, passes an object aperture at distance s Q , and then is focused onto 
photon-counting detector D\ by a collection lens. In arm-2 the idler propagates freely over a 
distance d2 from the output plane of the source to a point- like photon counting detector D2 ■ 



(I) Arm-1 (source to object): 

The optical transfer function or Green's function in arm-1, which propagates the field from 
the source plane to the object plane, is given by: 

g(n s ,u Js ;p 1 ,z 1 = di + -so) - e^ 21 / dp t [ dp s { Z^L e <«W. e *£&\?.-* I 2 } 

J lens J source ZTTCai > 

x e -<^r|ftl a f Z^le^l^-^I'V (22) 
I 2ircs J 

where p s and pi are the transverse vectors defined, respectively, on the output plane of the source 
and on the plane of the imaging lens. The terms in the first and second curly brackets in Eq. (22) 
describe free space propagation from the output plane of the source to the imaging lens and from 
the imaging lens to the object plane, respectively. Again, e* 2cd i ' ^ s ~^ ' and e l2c3 ° ' P' - ? 1 1 are the 
Fresnel phases as defined in the Appendix. Here the imaging lens is treated as a thin-lens, and the 
transformation function of the imaging lens is approximated as a Gaussian, Z(|pz|, /) = e l5 =? ' Pl ' . 

(II) Arm-2 (from source to image) : 

In arm-2, the idler propagates freely from the source to the plane of D2, which is also the plane 
of the image. The Green's function is 

fffa, Wii p 2 , z 2 = d 2 ) = e 1 ^ f dp' s e *i=fcl I 2 e ««-£ (23) 

Z7TCU2 J source 

where p' s and P2 are the transverse vectors defined, respectively, on the output plane of the source 
and the plane of photodetector D 2 ■ 
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(Ill) ^(pi,p 2 ) and G^ 2 \pi, p 2 ) (object plane - image plane): 

For simplicity, in the following calculation we consider degenerate (uj s — loi — to) and collincar 
SPDC. The effective transverse biphoton wavefunction ^(pi,p 2 ) is then evaluated by substituting 



the Green's functions g(H s , to; p\, Z\) and to; p 2 , z 2 ) into Eq. (21), 
^(pi,p 2 ) oc J dn s dKiS(U s + Hi) g(H s , u; pi, z{) g(Hi, w; p 2 , z 2 ) 

Ke <f (».+«*) [ dK a dKi6(K a + Ki) { dpi { dp s e^'^e^ 1 



1 u.< 



x e -*£f\p<\ j£k\pi-fc\ / dp' s e^' p = e* 5 * 1 "* - * 1 (24) 



where all the proportionality constants have been ignored. After completing the double integral of 
cK, and dHj 

dH s dKi 5(H S + Hi) e ilC '^' e iR '- p '° ~ 5(p s - fy, 



Eq. (24 1 becomes 

./ Zens ./ source 

Next, we complete the integral for dp s , 

V(pi,p 2 ) oc e ^ (so+5l) / dp /e ^ |pW!l V l w^lV^ |pWl12 , (25) 



where we have replaced d\ + d 2 with Si (as depicted in Fig. [8]). Although the signal and idler prop- 
agate in different directions along two optical arms, interestingly, the Green function in Eq. (25) 
is equivalent to that of a classical imaging setup, as if the field is originated from a point p\ on 
the object plane and propagated the lens and then arrived at point p 2 on the imaging plane. The 
mathematics is consistent with our previous qualitative analysis of the experiment. 

The finite integral on dpi yields a point-to- "spot" relationship between the object plane and 
the image plane that is defined by the Gaussian thin-lens equation 



./>..,. \s n c m / 



R w 

lens ^ s o C 

If the integral is taken to infinity, by imposing the condition of the Gaussian thin-lens equation 
the effective transverse biphoton wavefunction can be approximated as a S function 

*(pi, P2) ~ 5(Pi + fe/m) ~ 8{p + pi/m), (27) 

where we have replaced p\ and p 2 with p and pi, respectively, to emphasize the point-to-point 
EPR correlation between the object and image planes. To avoid confusion with the "idler" we 
have used pi to label the image plane. 

We now include an object-aperture function, a collection lens and a photon counting detector 
Di into the optical transfer function of arm-1 as shown in Fig. [2] The collection- lens— D\ package 
can be simply treated as a "bucket" detector. The "bucket" detector integrates the biphoton 
amplitudes ^{p ,p 2 ), which are modulated by the object aperture function A(p ) into a joint 
photodetection event. This process is equivalent to the following convolution 



#1,2 oc / dp \A(p )\ 2 \m(p ,p 2 )\ 

J object 



\A{p 2 /m)\ 2 = \A(pj/m)\ 2 . (28) 
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Again, D2 is scanned in the image plane (p2 = pi)- A ghost image of the object is thus reproduced 
on the image plane by means of the joint-detection between the point-like-detector D2 and the 
bucket detector D\. 

The physical process corresponding to the above convolution is rather simple. Suppose the 
point detector D2 is triggered by an idler photon at a transverse position of p\ in a joint-detection 
event with the bucket detector Di which is triggered by the signal twin that is either transmitted 
or reflected from a unique point p on the object plane. This unique point-to-point determination 
comes from the non-factorizable correlation function 8{p Q + pi/m). Now, we move D2 to another 
transverse position p' I and register a joint-detection event. The signal photon that triggers Di 
must be either transmitted or reflected from another unique point p' a on the object plane which 
is determined by S(p' Q + p'i/m). The chances of receiving a joint detection event at p\ and at 
p'l would be modulated by the values of the aperture function A{p ) and A{p' Q ), respectively. 
Accumulating a large number of joint-detection events at each transverse coordinates on the image 
plane, the aperture function A(p a ) is thus reproduced in the joint-detection as a function of p\. 



The observation of type-one ghost imaging has demonstrated a non-factorizable point-to-point 
EPR correlation between the object and image planes. This point-to-point correlation is the result 
of a constructive-destructive interference between biphoton amplitudes, 

G (2) (p ,p/) = \ dK s dn i b{K s + >i i )g{K s ,p )g{ii i ,pi) = somb 2 ( — —\p + —\) . (29) 



In this view we consider the ghost imaging experiment of Pittman et al. a realization of the 1935 
EPR gedankenexperiment |21j |22j . 

Classical theory has difficulties when facing type-one ghost imaging phenomenon. In the 
classical theory of light, a joint measurement between two photodetectors D\ and D2 measures the 
statistical correlation of intensity fluctuations, 

G (2) (pi,P2) - (I(pi)I(p 2 )) = I(pi)I(p2) + (AI(p 1 )AI(p 2 )). (30) 
Therefore, the point-to-point image-forming correlation is considered as a result of the statistical 



correlation of intensity fluctuations between the object and the image planes. Comparing Eq. (30 1 



which has a constant background I{pi)I{p2), with Eq. (29 1, which has a zero background, the 
mean intensities I(p\) and I(p2) must be zero, otherwise the result would lead to non-physical 
conclusions. The measurements, however, never yield zero mean values of I(p\) and I(p2) under 
any circumstances. In fact, the individual-detector counting rates of D\ and D2 were monitored 
in the experiment of Pittman et al. with much greater value than that of the coincidences. It is 
clear that the classical theory of statistical correlation of intensity fluctuations does not reflect the 
correct physics behind type-one ghost imaging. 



4 Type-two ghost imaging with chaotic radiation 

In this section we discuss the physics of type-two ghost imaging. The near-field lensless ghost 
imaging with chaotic radiation was first demonstrated by Scarcelli et al. in the years from 2005 to 
2006 [5 j [6 following their experimental demonstration of two-photon interference of chaotic light 
[7J. The schematic experimental setup of their 2006 demonstration is shown in Fig. [9j Radiation 
with a narrow spectral bandwidth Aw of a few millimeters diameter from a chaotic pseudothcrmal 
source [23] was equally divided into two by a 50% — 50% non-polarizing beam-splitter. In the 
reflected arm, a double-slit with slit separation b = 1.5 mm and slit width a — 0.2 mm, was 
placed at a distance z 2 = 139 mm from the source and a bucket detector D 2 was placed just 
behind the object. In the transmitted arm a point detector D\ was scanned in the transverse 
plane of z\ — Z2- Scarcelli et al tested two different joint detection schemes, namely the photon 
counting coincidence circuit and the standard HBT correlator. In the photon counting regime two 
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Figure 9: Near- field lcnslcss ghost imaging with chaotic light demonstrated in 2006 by Scarcelli 
et al. Di is a point-like photodetector that is scannable along the xi-axis. The joint-detection 
between D\ and the bucket detector D 2 is realized either by a photon-counting coincidence 
counter or by a standard HBT linear multiplier (RF mixer). In this measurement D 2 is fixed 
in the focal point of a convex lens, playing the role of a bucket detector. The counting rate or 
the photocurrent of Di and D 2l respectively, are measured to be constants. Surprisingly, an 
image of the 1-D object is observed in the joint-detection between D\ and D 2 by scanning D\ 
in the plane of z\ = z 2 along the cci-axis. The image, is blurred out when z\ ^ z 2 . There is 
no doubt that thermal radiations propagate to any transverse plane in a random and chaotic 
manner. There is no lens applied to force the thermal radiation "collapsing" to a point or 
speckle cither. What is the physical cause of the point-to-point image-forming correlation in 
coincidences? 



Geiger mode avalanche photodiodes were employed for single-photon detection. In the bright light 
condition, two silicon PIN diodes were used with a standard analog HBT linear multiplier. The 
bucket detector D 2 was simulated by using a short focal length lens (/ = 25mm) to focus the 
light coming from the object onto the active area of the detector while the point detector D 2 was 
simulated by a pinhole like aperture. After a large number of reaped measurements for different 
experimental schemes and conditions, Scarcelli et al reported the following observations. 



Observation (1): A typical measured ghost image of the double-slit is shown in Fig. 10 The 
measured curve reports the joint-detection counting rate between Di and D 2l or the output current 
of a HBT linear multiplier, as a function of the transverse position of the point detector Dx along 
X\ axis. Notice, in Fig. [10] the constant background has been removed from the correlation. 

Observation (2): The measured contrasts vary significantly under different experimental schemes 



Object Image 

III 

-1.5 -1 -0.5 0.5 1 1.5 
Position of D\ (mm) 

Figure 10: The double-slit and its ghost image. Notice, the constant background has been 
removed from the correlation. 
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and conditions. It was found that the image contrast can achieve ~ 50% in photon counting 
measurement if no more than one joint-detection event occurring within the time window of the 
coincidence circuit. 50% is the maximum image contrast we expect for thermal light ghost imaging. 

Observation (3): To achieve less than one joint-detection event per coincidence time window, 
weak light source is not a necessary condition. It can be easily achieved under bright light condition 
by using adjustable ND-filters with Di and D 2 - 

l\ t Secondary 




17.2 mm 



Figure 1 1 : Unfolded schematic experimental setup of a secondary image measurement of the 
ghost image and the measured secondary images. By using a convex lens of focal length /, 
the ghost image is imaged onto a secondary image plane, which is defined by the Gaussian 
thin-lens equation, l/s a + 1/sj = 1//, with magnification m = —Si/s Q . This setup is useful 
for distant large scale ghost imaging applications. 

To confirm the observations are imaged images, and not "projection shadows" , Scarcelli et al. 
made two additional measurements. In the first measurement, photodetector D\ was moved away 
from the ghost image plane of Z\ — z%. Whether moved in the direction of Z\ > z 2 or z\ < z%, the 
ghost image became "blurred" . The measurement also showed that the depth of the image is a 
function of the angular size of the thermal source: a larger angular sized source (opening angle AO 
relative to the photodetectors) produces sharper image with shorter image depth. In the second 
measurement, Scarcelli et al. constructed a secondary imaging system, illustrated schematically in 
Fig. 1 1 1 1 By using a convex lens of focus length / the ghost image is imaged onto a secondary image 
plane, which is defined by the Gaussian thin-lens equation, l/s a + = 1//, with magnification 
m = — Sj/s D . In this measurement the scanning photodetector Di is placed on the secondary 
imaging plane. The secondary image of the ghost image is observed in the joint-detection between 
D\ and D 2 by means of either a photon-counting coincidence counter or a HBT linear multiplier. 

4.1 What is the physical cause of chaotic light ghost imaging? 

It is the partial point-to-point correlation between the object plane and the image plane that 
makes ghost imaging with thermal light possible. Similar but different from classical imaging and 
type-one ghost imaging, mathematically, type-two ghost imaging is the result of a convolution 
between the aperture function | ^4 ( />2 ) | 2 and a (^-function like partial point-to- "spot" correlation 
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function 

i?i2 0c/ dp 2 \A(p 2 )\ 2 [1 + somb 2 ( ^ - % (31) 

t' object 

in 2-D, where A# is the angular diameter of the radiation source viewed from the photodetector, 
Pi and p*2 are the transverse coordinates on the object plane and the image plane, respectively, or 

i?i 2 oc f dx 2 \A(x 2 )\ 2 [1 + sinc 2 ( — ^ - X2) )] . (32) 



object 



in 1-D. For a chosen wavelength, the spatial resolution of the ghost image is determined by the 
angular diameter of the light source: the larger the size of the source in transverse dimensions, 
the higher the spatial resolution of the lensless ghost image. The point-to- "spot" image- forming 
functions in Eqs. (31 I and (32 1 have been verified experimentally by Scarcclli et al. 

The physical process corresponding to the convolution of Eq. (31 1 and (32 1 is similar to that 
of the type-one ghost imaging. Suppose the point detector D\ or a CCD clement is triggered by a 
photon at a transverse position of p\ in a joint-detection event with the bucket detector D 2 which 
is triggered by another photon that is either transmitted or reflected from the object. According 
to Eq. (31) and (32), under condition of Z\ = z 2 , the photon from the object would have twice 
greater chance to be found at p bj = pi- Now, we move D\ to another transverse position p\, or 
locate another CCD element at p\ for joint-detection. The photon that triggers D 2 would have 
twice greater chance of been located at p' bj = P'i- The probabilities of receiving a joint detection 
event at p\ = /? bj and at p\ = p' bj would be modulated by the values of the aperture function 
A(pobj) a- n d A(p' ob A, respectively. Accumulating a large number of joint-detection events for each 
transverse coordinates pi, or for each CCD element in the image plane, a 50% contrast aperture 
function A(pi) = A{p a \ s] ) is thus reproduced in the joint-detection as a function of Pij^] 

To achieve thermal light ghost image with 50% contrast, we need (1) randomly distributed 
radiations on the object plane and on the image plane, respectively; and (2) for any photoelectron 
event at p\ there exists a unique corresponding point p bj = p\ on the object plane which has twice 
chance of observing another photoelectron event jointly and simultaneously. There is no doubt that 
random and chaotic radiation would propagate to any transverse plane in a random and chaotic 
manner. Therefore, condition (1) is satisfied automatically for chaotic thermal radiation. However, 
it is not easy to understand condition (2). We have been asking ourself a question since the first 
observation of lensless thermal light ghost image: what is the physical cause of the non-factorizable 
partial point-to-point image- forming function of 1 + S( p\ — /9 bj )? There seems no reason to have 
such a statistical correlation for thermal light. Figure [l2| schematically illustrates this situation. 
To simplify the picture we assume the source in 1-D with a large number of independent point 
sub-sources randomly distributed from —6/2 to 6/2. Each point sub-source, such as the jth and the 
fcth sub-source, randomly radiates independent spherical waves to the object and image planes, 
respectively. Due to the chaotic nature of the source, these independent and incoherent sub- 
intensities simply add together yielding a constant total intensity spatially and temporally on any 
transverse plane. The more chaotic sub-fields that contribute to the intensity sum, the less value 
of AI/I is expected. For any two transverse planes, such as the object and the image planes in 
Fig. [9] each with independent and randomly distributed intensities, statistically, there is no reason 
to expect any spatial or temporal correlations. What is the physical cause that forces a twice large 
probability for the thermal radiation to jointly appear at p\ — p bj? 

In fact, we have been facing this question since 1956, after the discovery of Hanbury Brown 
and Twiss (HBT). The lensless ghost imaging setup looks similar to that of the historical HBT 
spatial interferometer which was used for measuring the angular size of distant stars. A significant 
difference is that the lensless ghost imaging measurement is in near-fielcQfor imaging purposes [5]. 

6 To observe thermal light ghost image with maximum 50% contrast requires achieving a necessary experimental 
condition: no more than one joint detection event within the coincidence time window. 

7 The concept of "near-field" was defined by Fresnel to be distinct from the Fraunhofer far-field. The Fresncl 
near-field is different from the "near-surface-field" which considers a distance of a few wavelengths from a surface. 
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Figure 12: A large number of independent point sub-sources, such as the jth and kth, are 
randomly distributed on the plane of a thermal source. These point sub-sources randomly 
radiate independent spherical waves to the object and image planes, respectively. Due to the 
chaotic nature of the source, these independent sub-intensities simply added together yielding 
a constant total intensity in space and in time on any transverse planes. 



The HBT experiment created quite a surprise in the physics community with an enduring debate 
about the classical or quantum nature of the phenomenon [TT][T2]. Figure 13 is a schematic of the 
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Figure 13: Schematic of the historical Hanbury Brown and Twiss experiment which measures 
the transverse spatial correlation of far-field thermal radiation. 

historical HBT experiment which measures the transverse spatial correlation of far-field thermal 
radiation. Performing the measurement in 1-D by scanning photodetectors D\ and/or D<i along 
the axes x\ and X2, the second-order transverse spatial correlation function (x\ , X2) was found 
to be 



1 



■nA6(x\ — 22) 
A 



}• 



(33) 



where A9 is the angular size of the star, A the wavelength of the radiation. The far-field HBT 



correlation of Eq. ( 33 1 has been interpreted as the result of classical statistical correlation of the 



intensity fluctuations 



I hh ) = ((h + &h){h + AZ 2 )) = hh + ( A/xA/a ) 



where 1\ and I2 are the mean intensities of the radiation measured by photodetectors D\ and D2 , 



respectively. The second term in Eq. (33 1, Iq sine [t:A6(xi — x 2 )/X], is phenomenologically inter- 



preted as the intensity fluctuation correlation (A/1A/2) in classical theory. For visible wavelengths 
and large values of A9 this function quickly drops from its maximum to minimum when x\ — a; 2 
moves from zero to a value such that AB{x\ — x%)l\ = 1. In this situation we effectively have a 
"point-to-point" relationship between the x\ and X2 axes: for each point on the x\ there exists 
only one point on the x 2 that may have a nonzero intensity fluctuation correlation. 
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Figure 14: A phenomenological interpretation of the historical HBT experiment. Upper: 
the two photodetectors receive identical modes of the far-field radiation and thus experience 
identical intensity fluctuations. The joint measurement of D\ and D% gives a maximum value 
of (Aii A/2) ■ Lower: the two photodetectors receive different modes of the far-field radiation. 
In this case the joint measurement gives (A/1A/2) = 0. Unfortunately, this hand-waving 
interpretation does not reflect the correct physics in the case of A9 7^ 0. For a finite angular 
sized source, there is no chance, at least realistically, for Di and D 2 to receive radiation from 
a single radiation mode only. Nevertheless, the above theory has convinced us to believe that 
the observation of the intensity fluctuation correlation only takes place in the far-field zone of 
the thermal source. 
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The well-accepted interpretation of the HBT phenomenon is the following: in HBT the mea- 
surement is taken in the far- field zone of the radiation source, which is equivalent to the Fourier 
transform plane. When D\ (D 2 ) is scanned in the neighborhood of x\ = x 2 , the two detectors 
measure the same mode of the radiation field. The measured intensities have the same fluctuations 
and yield a maximum value of (AI1AI2) ■ The two upper curves of I(t) in Fig. 15 schematically 
illustrate this situation. When the two photodetectors move apart from X\ = x 2 , D\ and D 2 
measure different modes of the radiation field. In this case, the measured two modes may have 
completely different fluctuations. The measurement yields (A/1A/2) = and gives (I1I2) = Iih- 



This situation is illustrated in the two lower curves of I(t) in Fig. 15 Unfortunately, this hand- 
waving interpretation does not reflect the correct physics in the case of AO ^ 0. For a finite 
angular sized source, there is no chance, at least realistically, for Di and D 2 to receive radiation 
from a single radiation mode only. Nevertheless, the above theory has convinced us to believe that 
the observation of the intensity fluctuation correlation only takes place in the far-field zone of the 
thermal source. 




Figure 16: Modified near-field HBT measurement - an unfolded Klyshko picture of the setup. 
Assuming a large sized disk-like near-field chaotic source, each point on the disk can be 
considered as an independent sub-source. It is easy to see that (1) D\ and D 2 are capable of 
receiving radiation from a large number of sub-sources; and (2) D\ and D 2 have more chances 
to be triggered jointly by radiations from different sub-sources; (3) The ratio between the 
joint-detections triggered by radiation from a single sub-source and from different sub-sources 
is roughly N/N 2 = 1/N in any transverse position of D 1 and D 2 . 



What will happen if we move the photodetectors Di and D 2 to the "near-field" as shown in 
the unfolded schematic of Fig. \lEf Does this hand-waving argument still predict the point-to- 
point correlation in this situation? We consider a disk-like thermal source with a large number 
of independent and randomly radiating point sub-sources and assume the radiations coming from 
the same sub-source have the same intensity fluctuation, and the radiations coming from different 
sub-sources have different intensity fluctuations. It is easy to see that in the near-field, (1) each 
photodetector, D\ and D 2 , is capable of receiving radiations from a large number of sub-sources; 
and (2) D\ and D 2 , have more chances to be triggered jointly by radiation from different sub- 
sources; (3) The ratio between the joint-detections triggered by radiation from a single sub-source 
and from different sub-sources is roughly N/N 2 — 1/N in any transverse position of D\ and D 2 . 
For a large value of N the contribution of joint-detections triggered by radiation from a single sub- 
source in any transverse position of D\ and D 2 has the same negligible value (AIiAI 2 ) /Iil 2 ~ 0. 
Following the above philosophy, the near-field G^ 2 \p\,p 2 ) should be a constant for any chosen 
transverse coordinates p\ and p 2 . The experimental observations, however, have shown a different 
story. 

The nontrivial near-field point-to-point correlation was experimentally observed in a modified 
HBT experiment by Scarcelli et al. in 2005 before the near-field lensless ghost imaging demon- 
stration. The modified HBT has a similar experimental setup as that of the historical HBT of 
Fig. 13 except replacing the distant star with a near-field disk-like chaotic source. This light 
source has a considerably large angular diameter from the view of the photodetectors Di and D 2 . 
The point photodetectors Di and D 2 are scannable along the axes of X\ and x 2 , respectively. The 
frequency bandwidth Auj of this thermal source is chos en to be narrow enough to achieve ^ps 
correlation width of G^ 2 ^(ti — t 2 ) which is shown in Fig. 17 This means to change G^ 2 ) from its 
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Figure 17: G^(h - t 2 ) of a chaotic source. The temporal correlation width is measured 
~ 0.5/xs, which means that unless t\ — t% > 0.5/xs, the value of G^ 2 '{t\ — t 2 ) will stay at the 
neighborhood of its maximum. 

maximum (minimum) value to minimum (maximum) value requires a few hundred meters optical 
delay in the arm of either Di or D 2 . The transverse intensity distributions were examined before 
the measurement of transverse correlation. The counting rate (weak light condition) or the output 
photocurrent (bright light condition) of each individual photodetector was found to be constant, 
i.e., I{p\) ~ constant and I(p2) ~ constant by scanning D\ and D 2 in the transverse planes of 
Z\ — zq and z 2 = Zo, respectively. There is no surprise to have constant I(pi) and I{p 2 )- The 
physics has been clearly illustrated in Fig. [12] By using this kind of chaotic source, Scarcelli et al. 
measured the 1-D near-field normalized transverse spatial correlation of g^ 2 \x\ — x 2 ) by scanning 
D\ in the neighborhood of x% — x 2 . The measurements confirmed the point-to- "spot" correlation 
of 



g( 2 \xi — x 2 ) ~ 1 + sine 2 [ 



nA9(xi — x 2 ) 
A 



(34) 



where, again, AO is the angular diameter of the near-field disk-like chaotic source. It is worth 
emphasizing that g( 2 \x\ — x 2 ) dependents on x\ — x 2 only. Taking x\ — x 2 = constant, g^ 2 \x\ — x 2 ) 
is invariant under the displacements of transverse coordinates. 

A simplified summary of the experimental observation is shown in Fig. 



18 (1) In the upper 



figure, D\ and D 2 are placed at equal distances from the source and aligned symmetrically on the 
optical axis. The normalized joint-detection, or the value of g^ achieved its maximum of ~ 2. 
(2) In the middle figure, D\ is moved up a few millimeters to a non-symmetrical position, the 
normalized joint-detection, or the value of g^ is measured to be ~ 1. (3) In the lower figure, 
D 2 is moved a few millimeters up to a symmetrical position with respect to D\. The normalized 
joint-detection, or the value of g^ turned back to its maximum of ~ 2 again. 

It is easy to see that the classical theory of statistical correlation of intensity fluctuations is 
facing difficulties in explaining the experimental results. In near-field D\ and D 2 receive the same 
large number of modes at any p\ and p 2 . In the spirit of the traditional interpretation of HBT, 
there seems no reason to have a different i nte nsity fluctuation correlation between pi = p 2 and 
Pi 7^ Pi f° r the G^ 2 ) function shown in Fig. 
maximum value of g^ 
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In the upper measurement, we have obtained the 
2 at Z\ — z 2 and x\ — x 2 , which indicates the achievement of a maximum 
intensity fluctuation correlation as shown in Fig. 17 with \t\ — 1 2 \ ~ 0. In the middle measurement, 
g( 2 ) r^i 1 indicates a minimum intensity fluctuation correlation by moving D\ a few millimeters up, 
which means the intensities measured by D\ and D 2 must have different fluctuations. In the lower 
measurement D 2 is moved up a few millimeters to a new symmetrical position with respect to Di, 
the measurements obtain g^> ~ 2 again. The intensities measured by D\ and D 2 must have same 
fluctuations again. What is the physical cause of the changes of the intensity fluctuations then? 
Remember the G^'{t\ — t 2 ) function has a width of ~ 0.5/z. 

For half a century since 1956, it has been believed that the HBT correlation is observable in 
the far-field only. It was quite a surprise that in 2005 Scarcelli et al. successfully demonstrated 
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Figure 18: Schematic of the near-field spatial correlation measurement of Scarcelli et al. 
Upper: ~ 2, where D\ and D2 are placed at equal distances from the source and aligned 
symmetrically on the optical axis. In the spirit of the traditional interpretation of HBT 
the intensities measured by D\ and D2 must have same fluctuations as shown in the figure. 
Middle: g^ ~ 1, where D± is moved up a few millimeters to an asymmetrical position. In 
the spirit of the traditional interpretation of HBT the intensities measured by D\ and D% 
must have different fluctuations. Lower: g^ ~ 2, where D% is moved up to a symmetrical 
position with respect to Di, again. In the spirit of the traditional interpretation of HBT the 
intensities measured by D\ and D2 must have same fluctuations again. What is the physical 
cause of the changes of the intensity fluctuations then? Remember the G^ 2 '(ti — £2) function 
has a width of ~ 0.5/i, see Fig. [TTj 

a near-field point-to-point transverse correlation of chaotic light, indicating that the nontrivial 
HBT spatial correlation is observable in the near-field and is useful for reproducing ghost images 
in a nonlocal manner]^] The experiment of Scarcelli et al. raised a question: "Can two-photon 
correlation of chaotic light be considered as correlation of intensity fluctuations?" [5] At least, this 
experiment suggested we reexamine the relationship between the quantum mechanical concept of 
joint-detection probability with the classical concept of intensity fluctuation correlation. It seems 
that jointly observing a pair of photons at space-time point (ri,ti) and (rjjta) is perhaps only 
phenomenologically connected but not physically caused by the classical statistical correlation of 
intensity fluctuations. The point-to-point image-forming correlation is more likely the result of 
an interference. In the view of two-photon interference, far-field is not a necessary condition for 
observing the partial point-to-point correlation of thermal light. Furthermore, it is quite common 
in two-photon interference type experiments to observe constant counting rates or intensities in 
individual photodetectors D\ and D2, respectively, and simultaneously observe nontrivial space- 
time correlation in the joint-detection between D\ and D2. These observations are consistent with 
the quantum theory of two-photon intcrfcrometry [5] . 

4.2 Quantum theory of thermal light ghost imaging 

According to the quantum theory of light, the observed partial point-to-point image- forming 
correlation is the result of multi-photon interference. In Glauber's theory of photo-detection |20j . 
an idealized point photodetector measures the probability of observing a photo-detection event at 

8 We cannot help but stop to ask: What has been preventing this simple move from far-field to near-field for half 
a century? 
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space-time point (r, t) 

G (1) (r,i) = tr{pE^(r,t)E^{r,t)}, (35) 

where p is the density operator which characterizes the state of the quantized electromagnetic 
field, E^'(r,t) and i?' + '(r,t) the negative and positive field operators at space-time coordinate 
(r,i). The counting rate of a point photon counting detector, or the output current of a point 
analog photodetector, is proportional to G^'(r,t). A joint-detection of two independent point 
photodetectors measures the probability of observing a joint-detection event of two photons at 
space-time points (ri,ti) and (r 2 ,t 2 ) 

G (2 \r 1 ,t 1 ;r 2 ,t 2 )=tr{pE^\r 1 ,t 1 )E^\r 2 ,t 2 )E^\v 2 ,t 2 )E^\r 1 ,t 1 )}, (36) 

where (rj,tj), j = 1,2, is the space-time coordinate of the jth photo-detection event. The coinci- 
dence counting rate of two photon counting detectors, or the output reading of a linear multiplier 
(RF mixer) between two photodetectors, is proportional to G^- 2 '(ri,ti;r 2 ,t 2 ). To calculate the 
partial point-to-point correlation between the object plane and the image plane, we need (1) to 
estimate the state, or the density matrix, of the thermal radiation; and (2) to propagate the field 
operators from the radiation source to the object and the image planes. We will first calculate the 
state of thermal radiation at the single-photon level for photon counting measurements to explore 
the physics behind ghost imaging as two-photon interference and then generalize the result to any 
intensity of thermal radiation. 

We assume a large transverse sized chaotic source consisting of a large number of independent 
and randomly radiating point sub-sources. Each point sub-source may also consist of a large 
number of independent atoms that are ready for two-level atomic transitions in a random manner. 
Most of the time, the atoms are in their ground state. There is, however, a small chance for each 
atom to be excited to a higher energy level E 2 (AE 2 ^ 0) and later return back to its ground 
state Ei. It is reasonable to assume that each atomic transition generates a field in the following 
single-photon state 

|*)~|0) + eX;/(k J *)4 > J0>. ( 37 ) 

where |e| <C 1 is the probability amplitude for the atomic transition, /(k, s) = ("fk.s |^) is the 
probability amplitude for the radiation field to be in the single-photon state of wave number k and 
polarization s: | ^k.s ) = | lk.s ) = s | ). For this simplified two-level system, the density matrix 
that characterizes the state of the radiation field excited by a large number of possible atomic 
transitions is thus 



p = II { 1 ) + e E z( k > s ) e ~ lu)t03 «L 1 ) } ( 38 ) 

toj k,s 

x II {<°l + e *E /(k', s V w% M0K'„ 5 <} 

*ofc k',s' 

^{l») + CE /( k < s ) e "" t0j C I >] + e 2 [...} } 

toj k,s 

x {(01 +^1EE /(k', S ')e^' tofc (0|a k 'v] + e * 2 [...]}, 

tok k',s' 

where e~ luJtoj is a random phase factor associated with the jth atomic transition. Since |e| <C 1, it 



is a good approximation to keep the necessary lower-order terms of e in Eq. ( 38 1 . After summing 
over toj (tok) by taking into account all its possible values we obtain 

/ 5^|0)(0| + | e | 2 ^|/(k, S )| 2 |l k , s )(l k , s | 

k,s 

+ l £ | 4 E E l/( k ' s )! 2 l/( k V)| 2 1 Ik.slk'.yXlk.slk',*' I + (39) 

k.s 
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Similar to our earlier discussion we will focus our calculation on the transverse correlation by 



assuming a narrow enough frequency bandwidth in Eq. ( 39 1 . In the experiments of Scarcelli et 



al. the coherence time of the radiation was chosen ~/xs, the maximum achievable optical path 
differences ^ps by the scanning of D\ and D 2 , and the response time of the photodetectors is 
much less than the coherence time. The transverse spatial correlation measurement is under the 
condition of achieving a maximum temporal coherence of G^ 2 \t\ — t 2 ) ~ 2 during the scanning 
of D\ and D 2 at any p\ and p 2 . In the photon counting regime, under the above condition, it is 
reasonable to model the thermal light in the following mixed state 



|o>(o| + | e | 2 ^a+(^)|o>(o|«(^) + | e | 4 ^^at(^)at(«') |o>(o|«(^')a(«)- 



(40) 



Basically we are modeling the light source as an incoherent statistical mixture of single-photon 
states and two-photon states with equal probability of having any transverse momentum. The 
spatial part of the second-order coherence function is thus calculated as 

G (2) (Pi ,*i 5 P2 , z 2 ) = tr[ p E ( -*> (pi , Zl (p 2 , z 2 )EM (p 2 , z 2 )E^ (pi ,z x )\ 

= Y,{UU^-\p 1 ,z 1 )E^\p 2 ,z 2 )E^\^,z 2 )E^\p 1 ,z l )\l H l w ) 



^2 \^ft,K'(Pl,Zl,P2,Z 2 )\ 



(41) 



where we have defined an effective two-photon wavefunction in transverse spatial coordinates 

■*r,w{pi,zi\P2,z 2 ) = (0|£;W(P2,Z2)^ (+) (pi^i)|1k1k')- (42) 
The transverse part of the electric field operator can be written as 

E {+) (pj ,Zj)oc^2 9j (Pj , zj ; k) a(it) , (43) 



again, gj(pj,Zj\ it) is the Green's function. Substituting the field operators into Eq. (42 1 we have 
^K,K'{PuZi;p 2 ,z 2 ) = ^[g 2 {p2, z 2 \K)g 1 {p 1 , zx\k') + g 2 {p^, z 2 \K')gi{pi, zr,H)] (44) 



and 



G ( &(p 1 ,zi;p 2 ,z 2 ) = ^ I -j=[g 2 (p 2 , z 2 ; Z)gi(p\, Z\\ k') + g 2 (p 2 , z 2 ; tf)g x (p u z x ; k)] 



(45) 



representing the key result for our understanding of the phenomenon. Eqs. (44 1 and (45) indicates 



an interference between two alternatives, different yet indistinguishable, which leads to a joint 
photo-detection event. This interference phenomenon is not, as in classical optics, due to the 
superposition of electromagnetic fields at a local point of space-time. This interference is the 
result of the superposition between g 2 (p 2 , z 2 ; K)gi(pi, Z\\ k!) and g 2 {p 2 , z 2 ; K,')gi(pi, Z\\ it), the so- 
called two-photon amplitudes, non-classical entities that involve both arms of the optical setup 
as well as two distant photo-detection events at (pi,Z\) and (p 2 ,z 2 ), respectively. Examining the 
effective wavefunction of Eq. ( 44 ) , we find this symmitrized effective wavefunction plays the same 
role as that of the symmitrized wavefunction of identical particles in quantum mechanics. This 
peculiar nonlocal superposition has no classical correspondence, and makes the type-two ghost 
image turbulence-free, i.e., any phase disturbance in the optical path has no influence on the ghost 
image [24]. Fig. 19 schematically illustrates the two alternatives for a pair of mode k and it to 
produce a joint photo-detection event: itl x it 2 and it2 x itl. The superposition of each pair 
of these amplitudes produces an individual sub-interference-pattern in the joint-detection space 
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Figure 19: Schematic illustration of two-photon interference: a superposition between two- 
photon amplitudes g2{p 2 , z 2^)gi{pi, z i] and <72(P2j z i\ k')9i{Pi, z \\ It is clear that the 
amplitudes (72 (f>2, z 2', k)9i(Pi, Z\\ i&) and 32 (p2, z 2', K r )gi(pi, z\ \ k) will experience equal optical 
path propagation and superpose constructively when D\ and D2 are located at pi ~ p2 and 
z\ ~ Z2. This nonlocal superposition has no classical correspondence. 



of (pi, iii Pi, ^2) • A large number of these sub-interference-patterns simply add together 
resulting in a nontrivial G^ 2 \p\, z\\ p 2 , z 2) function. It is easy to see that each pair of the two- 
photon amplitudes, illustrated in Fig. |19| will superpose constructively whenever D\ and D2 are 
placed in the positions satisfying pi ~ p*2 and z\ ~ z^\ and consequently, G^ 2 \p\, z\ \ P2, Z2) achieves 
its maximum value as the result of the sum of these individual constructive interferences. In other 
coordinates, however, the superposition of each individual pair of the two-photon amplitudes may 
yield different values between constructive maximum and destructive minimum due to unequal 
optical path propagation, resulting in an averaged sum. 

Before calculating G^ 2 \p\, Z\\p2, 22) we examine the single counting rate of the point photode- 
tectors Di and D2 which are placed at (pi,Zi) and (p 2 ,Z2), respectively. With reference to the 
experimental setup of Fig. [9j the Green's function of free-propagation is derived in the Appendix 



3i(Pi,zi;k) = 

92{P2-,Z2\K) 



dps I 
dp' s { 



-ioj e 



3 |2 - 



2ttc Z\ 



2irc Z2 



-IK-ps 



where p s is the transverse vector in the source plane, and the field has propagated from the source 
to the pi plane and p*2 plane in arms 1 and 2, respectively. The sing le d etector counting rate or 
the output photocurrent is proportional to G^(r,t) as shown in Eq. (35 1, 



G^{p v z 3 ) = ir{p^-)(^ l2i )^ (+) (Pi,^-)} 

oc 1 &{k) E(-\p , Zj (p} , Zj ) at (k) 1 ; 



constant, 



(46) 



where j — 1,2 indicating the jth photodetector. 

Although G^(p 1: zi) and G^{p 2 ,z 2 ) are both constants, G^ 2 \pi, Z\\ p 2 , Z2) turns to be a 
nontrivial function of (pi, z\) and (p 2 , Z2), 



G (2) (pi,zi;p2,z 2 ) = ^2 ^-^=[g 2 (p2,Z2;K)gi{pi,z 1 ;K') + g 2 (p 2 , z 2 ;Fi')gi(pi, zi;k)] 
= G (1) (pi^i)G (1) (p2,^)+ \G^\p 1) z 1 ;p2,Z2)\ 2 ) 



(47) 



where 



|G (1) (pi,zi;p 2 ,z 2 )| =| / dn gl{pi,zi;n) g 2 {p2,z 2 ;K) 

dp s e~^ Zl e" ! 



rlP2-AI 2 | 2 
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If we choose the distances from the source to the two detectors to be equal (z% = Z2 = d), the 
above integral of dp s yields a point-to-point correlation between the transverse planes z\ — d and 
z 2 = d, 



G$(pi;p 2 ) I 2 oc 



dp s e 



R u) 

= somb 2 [— — |pi - p 2 |] ~ S(pt - p 2 ), (48) 



where the 5-function is an approximation by assuming a large enough thermal source of angular 
size A9 ~ R/d and high enough frequency w, such as a visible light source. The nontrivial G^ 
function is therefore, 



G (2) (pi;p 2 ) ~i + <J(pi- P2). 

In the ghost imaging experiment, the joint-detection counting rate is thus 
i?i2<x / dp 2 \A{p 2 )\ 2 G {2 \p 1 -p 2 )~R () +\A{p 1 )\ 2 , 



(49) 



(50) 



where Rq is a constant and A(p 2 ) is the aperture function of the object. 

So far, we have successfully derived an analytical solution for ghost imaging with thermal 
radiation at the single-photon level. We have shown that the partial point-to-point correlation of 
thermal radiation is the result of a constructive-destructive interference caused by the superposition 
of two two- photon amplitudes, corresponding to two alternative ways for a pair of jointly measured 
photons to produce a joint-detection event. In fact the above analysis is not restricted to single- 
photon states. The partial point-to-point correlation of G^ 2 '(pi, p 2 ) is generally true for any order 
of quantized thermal radiation [25 . Now we generalize the calculation to an arbitrary quantized 
thermal field with occupation number from rtk.s = to r&k,s 3* 1 by keeping all higher order terms 
in Eq. ( 38 ) . After summing over toj and iofc the density matrix can be written as 



P=Y1 p M IWKWI' 



(51) 



where P{„} is the probability for the thermal field in the state 



IW) 



\n-k, a )\n v 



The summation of Eq. (51 1 includes all possible modes k, polarizations s, occupation numbers rik,. 



for the mode (k, s) and all possible combinations of occupation numbers for different modes in a 



set of {n}. Substituting the field operators and the density operator of Eq. (511 into Eq. ( |35| ) we 
obtain the constant G^(pj,Zj,tj), j — 1,2, which corresponds to the intensities I{p\, z\,ti) and 
I{P2,z 2 ,t 2 ), 

= ^2 / d K I dK'g*(p J1 z 1 ,tf,K)g J {p v z 1 ,t j -,K , )p {n} ({n}\ a(K)a\n') \{n}) 
{n} J J 



^n R n R , J dH \gj(pj,Zj,tf,Z) 



{«} 

constant 



(52) 



Although G^'(pi, zi,ti) and G^(p 2 , z 2 ,t 2 ) are both constants, substituting the field operators and 



the density operator of Eq. (51 1 into Eq. (36 1, we obtain a nontrivial point-to-point correlation 
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function of G^ 2 \p\] p^) at the two transverse planes z\ = d and z 2 = d, 
G {2) (pi;p2) 

= ^ J dti J dn' J dn" J diZ'" g{ (pi , z\ ; ti)g 2 {p 2 , z 2 ; K.')g 2 (p 2 , z 2 ; tf')gi(p\,zt',tf") 
W 

x p{„ } ({n}| a(K)a(K , )a t (K")a t (K'") \{n}) 

oc n g ng/ y / ^ J dtf' J dtf"gl(pi,zi',iz)gZ(f>2, z 2 ;K')g2{p2,Z2]tf')gi(pi,zi;K,''') 
M 

X P{n}(<W" + <W'£«'£"') 

= X! ?>{...«»...««/...} 



1 | 2 i 

^ [^(pl^K^O^,^;/*') + g 2 (p 2 , z 2 ; ^giipi, zi, k')] | | 



(x {l + somft 2 ^^^ -p 2 )]}. (53) 



It is clear that in Eq. ( |53[ ), the partial point-to-point correlation of thermal light is the result of a 
constructive-destructive interference between two quantum-mechanical amplitudes. We also note 



from Eq. (53 1 that the partial point-to-point correlation is independent of the occupation numbers, 
{n}, and the probability distribution, P{„j, of the quantized thermal radiation. 

It is interesting but not surprising to see that the effective two-photon wavefunction in bright 
light condition 

^k,k'(Pi,zi;p 2 , z 2 ) = [g 2 (p2, z 2 ; K)gi{pi, Zf, k') + g 2 (p 2 , z 2 ; tf)gi(p*u z^k)] 

is the same as that of weak light at single-photon level. In fact, the above effective wavefunction 
does play the same role in specifying two different yet indistinguishable alternatives for the two 
annihilated photons contributing to a joint-detection event of Di and D 2 , which implies that the 
partial point-to-point correlation is the result of two-photon interference in bright light condition. 
This nonlocal partial correlation indicates that a 50% contrast ghost image is observable at bright 
light condition provided registering no more than one coincidence event within the joint-detection 
time window. This requirement can be easily achieved by using adjustable ND-filters with D\ and 
D 2 . 

Quantum theory predicts and calculates the probability of observing a certain physical event. 
The output photocurrent of an idealized point photodetector is proportional to the probability 
of observing a photo-detection event at space-time point (r, t). The joint-detection between two 
idealized point photodetectors is proportional to the probability of observing a joint photo-detection 
event at space-time points (ri,ti) and {v 2 ,t 2 ). In most of the experimental situations, there exists 
more than one possible alternative ways to produce a photo-detection event, or a joint photo- 
detection event. These probability amplitudes, which arc defined as the single-photon amplitudes 
and the two-photon amplitudes, respectively, are superposed to contribute to the final measured 
probability, and consequently determine the probability of observing a photo-detection event or a 
joint photo-detection event. In the view of quantum theory, whenever the state of the quantum 
system and the alternative ways to produce a photo-detection event or a joint photo-detection 
event are determined, the result of a measurement is determined. We may consider this as a basic 
criterion of quantum measurement theory. 

4.3 A semiclassical model of nonlocal interference 

The multi-photon interference nature of type-two ghost imaging can be seen intuitively from 
the superposition of paired-sub-fields of chaotic radiation. Let us consider a similar experimental 
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setup as that of the modified HBT experiment of Scarcelli et al. . We assume a large angular sized 
disk-like chaotic source that contains a large number of randomly radiating independent point 
"sub-sources" , such as trillions of independent atomic transitions randomly distributed spatially 
and temporally. It should be emphasized that a large number of independent or incoherent sub- 
sources is the only requirement for type- two ghost imaging. What we need is an ensemble of 
point-sub-sources with random relative phases so that the sub-fields coming from these sub-sources 
are able to take all possible values of relative phases in their superposition. It is unnecessary to 
require the radiation source to have either nature or artificial intensity fluctuations at all. In this 
model, each point sub-source contributes to the measurement an independent spherical wave as 
a sub-field of complex amplitude Ej — aje lVj , where cij is the real and positive amplitude of the 
jth sub-field and tpj is a random phase associated with the jth sub-field. We have the following 
picture for the source: (1) a large number of independent point-sources distribute randomly on the 
transverse plane of the source (counted spatially); (2) each point-source contains a large number 
of independently and randomly radiating atoms (counted temporally); (3) a large number of sub- 
sources, either counted spatially or temporally, may contribute to each of the independent radiation 
mode (k, of) at Di and D 2 (counted by mode). The instantaneous intensity at space-time (rj-,i,-), 
measured by the jth idealized point photodetector Dj, j = 1,2, is calculated as 

I(v 3 ,t 3 ) = ETfattAEiTjttj) = J2 3f fo> *j) E E rn(rj,tj) 

I m 

= J2 Et^tfiE^tj) + £ *?fo>*i)£™fa>*i). ( 54 ) 

l—rn l^m 

where the sub-fields are identified by the index I and m originated from the / and m sub-sources. 
The first term is a constant representing the sum of the sub-intensities, where the Ith sub-intensity 
is originated from the ith sub-source. The second term adds the "cross" terms corresponding to 
different sub-sources. When taking into account all possible realizations of the fields, it is easy to 
find that the only surviving terms in the sum are these terms in which the field and its conjugate 



come from the same sub-source, i.e., the first term in Eq. ( 54 1 . The second term in Eq. ( 54 1 



vanishes if <pi — <p m takes all possible values. We may write Eq. (54 1 into the following form 

7(r J -,t J ) = (Z(r,t))+A/(r,t), (55) 

where 

(J(r,t)) = (£Et(r j ,t j )J2E Tn (r j ,t j ))=J2Et(r j ,t j )E l (r 3 :,t j ). (56) 

l in I 

The notation (...) denotes the mathematical expectation, when taking into account all possible 
realizations of the fields, i.e., taking into account all possible complex amplitudes for the large 
number of sub-fields in the superposition. In the probability theory, the expectation value of a 
measurement equals the mean value of an ensemble. In a real measurement, the superposition 
may not take all possible realizations of the fields and consequently the measured instantaneous 
intensity 7(r,i) may differ from its expectation value (I(r,t)} from time to time. The variation 
5I(r,t) turns to be a random function of time. The measured I(j,t) fluctuate randomly in the 
neighborhood of (I(r,t)} non-deterministically. 

In the classical limit, a large number of independent and randomly radiated sub-sources con- 
tribute to the instantaneous intensity I(rj,tj). These large number of independent randomly 
distributed sub-fields may have taken all possible realizations of their complex amplitudes in the 
superposition. In this case the sum of the cross terms vanishes, 

AI(r, t) = Y / m^t^E^ti) c 0, (57) 

therefore, 



Jfo.i,) ~ Y,Et(Ti,tj)Ei(T j ,t j ) - (J(r,i)>. 
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Now we calculate the second-order correlation function G( 2 )(ri,ii;r2, t 2 ), which is defined as 



GW(r 1 ,t 1 ;r 2 ,t 2 ) = < £,*(ri,ti) E k {r x ,t x )E^ (r 2 ,t 2 ) S m (r 2 ,t 2 ) ), (58) 

where the notation (...), again, denotes an expectation operation by taking into account all possible 
realizations of the fields, i.e., averaging all possible complex amplitudes for the sub-fields in the 
superposition. In the following calculation we only take into account the random phases of the sub- 
fields without considering the amplitude variations. Due to the chaotic nature of the independent 
sub-sources, after taking into account all possible realizations of the phases associated with the 
sub-fields, the only surviving terms in the summation are those with: (1) j = k,l = m, (2) 
7 = m, k = I. Therefore, G^fri 

j ti'i r 2, ^2) reduces to the sum of the following two groups: 

G< 2 >(ri,ii;r 2 ,t 2 ) - E* (r^h) E^r^t,) ]T Ef (r 2 , t 2 ) E^, t 2 ) 
3 1 

+ Y J E*{v u t 1 )E J {v 2l t 2 ) Y, £f(r 2 ,i2)£z(ri,ti)) 
3 1 

= <EE |^[^(ri,ti)^(r 2 ,t 2 )+^(ri,ti)^(r 2 ,i 2 )]| 2 ). (59) 



It is not difficult to see the nonlocal nature of the superposition shown in Eq. (|59j) . In Eq. ( 59 1 , 
G^(ri, ti;r 2 , t 2 ) is written as a superposition between the paired sub-fields Ej(r x , ti)_E/(r 2 , t 2 ) and 
Ei(ri,ti)Ej(r2,t2)- The first term in the superposition corresponds to the situation in which the 
field at Di was generated by the jth sub-source, and the field at D 2 was generated by the Ith 
sub-source. The second term in the superposition corresponds to a different yet indistinguishable 
situation in which the field at D\ was generated by the Ith sub-source, and the field at D 2 was 
generated by the jth sub-source. Therefore, an interference is concealed in the joint measurement of 
D\ and D 2 , which physically occurs at two space-time points (1*1, t\) and (r 2 ,i 2 ). The interference 



corresponds to \Ej\E%2 + Ei\E^\ . It is easy to see from Fig. 20 the amplitude pairs jl x 12 




Figure 20: Schematic illustration of J^j 1 \EjiE12 + EnEj2\ 2 ■ It is clear that the amplitude 
pairs jl x 12 with 11 x j'2, where j and I represent all point sub-sources, pair by pair, will 
experience equal optical path propagation and superpose constructively when D\ and D 2 are 
located at p\ — p 2 , Z\ ~ z 2 . This interference is similar to symmetrizing the wavefunction of 
identical particles in quantum mechanics. 

with 11 x j2, j'l x 1'2 with I'l x j'2, jl x 1'2 with I'l x j2, and j'l x 12 with 11 x j'2, etc., pair 
by pair, experience equal total optical path propagation, which involves two arms of D± and D 2 , 
and thus superpose constructively when D\ and _D 2 are placed in the neighborhood of p\ = p 2 , 
Z\ — z 2 . Consequently, the summation of these individual constructive interference terms will 
yield a maximum value. When p\ ^ p 2 , z\ — z 2 , however, each pair of the amplitudes may achieve 
different relative phase and contribute a different value to the summation, resulting in an averaged 
constant value. 
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It does not seem to make sense to claim a nonlocal interference between [(Ej goes to Di) 
x (Ei goes to D 2 )] and [(Ei goes to Di) x (Ej goes to D 2 )j in the framework of Maxwell's 
electromagnetic wave theory of light. This statement is more likely adapted from particle physics, 
similar to symmetrizing the wavefunction of identical particles, and is more suitable to describe the 
interference between quantum amplitudes: [(particle-j goes to D\) x (particle-Z goes to D 2 )] and 
[(particle-i goes to D{) x (particle-j goes to D 2 )}, rather than waves. Classical waves do not behave 
in such a way. In fact, in this model each sub-source corresponds to an independent spontaneous 
atomic transition in nature, and consequently corresponds to the creation of a photon. Therefore, 
the above superposition corresponds to the superposition between two indistinguishable two-photon 
amplitudes, and is thus called two-photon interference [9]. In Dirac's theory, this interference is 
the result of a measured pair of photons interfering with itself. 

In the following we attempt a near -fie ld calculation to derive the point-to-point correlation of 
G^ 2 \p\, z\ \ p 2 , z 2 ). We start from Eq. (|59| and concentrate to the transverse spatial correlation 



G (2) (p 1 ,z 1 ;p 2 ,z 2 ) = (J2J2\^l. E o^ z ^ E i^^ + E i ^i)Ej(p2,z 2 )]\ 2 }. (60) 
j 1 

In the near-field we apply the Fresnel approximation as usual to propagate the field from each sub- 
source to the photodetectors. G^ 2 '(pi, Z\\ p 2 , z 2 ) can be formally written in terms of the Green's 
function, 



ZV,P2,Z 2 ) 



dKdn\^j= [g(p 1 ,z 1 ,K)g(p 2l z 2 ,k') + g(p 2 , z 2 , K)g(pt,Zi,K')) 



(J dn \g(p 1 ,z 1 ,K)\ J dn' \g(p 2 , z 2 ,k')\ + \ / dn g* (pi, z x , k) g(p 2 , z 2 , k) \ ) (61) 
G«(ft, Zl )G^(p 2 , z 2 ) + \G^(p 1 ,z 1 ; p 2 ,z 2 )\ 



In Eq. (61 1 we have formally written G^ in terms of the first-order correlation functions G^\ 
but keep in mind that the first-order correlation function G^ and the second-order correlation 
function G^ represent different physics based on different measurements. Substituting the Green's 
function derived in the Appendix for free propagation 

g(p j ,Z j ,K) = — — ^— J dp a(po) e ^ Po > e 2 "o 

into Eq. (61 1, we obtain G^(pi, zi)G^(p 2 , z 2 ) ~ constant and 

|G«(p 1 ,z 1 ;/? 2 ,z 2 )| 2 cx |(— / dp Q a 2 (p )e-^^ e^^ 1 -^^^^ e l ^^-Po\ 2 }\ 2 . 



Z\Z 2 

Assuming a 2 (po) ~ constant, and taking z\ — z 2 = d, we obtain 

IgW^;^)! 2 oc I / dp a 2 (p )e-^^-^ 2 e *^alA-Al a 



oc \ e ~^(\Pi\ 2 -\P2\ 2 ) / dpva 2 (p )e l ^ { ^-P 2) -P°\ 2 



R LU 



cx somb 2 [ — - \px -p 2 \], (62) 
L d " J 



where we have assumed a disk-like light source with a finite radius of R. The transverse spatial 
correlation function G^(p\;p 2 ) is thus 

G (2) (|/?i ~M) = 7 o [! + somb2 ( f ~ c I* " M )] ■ (63) 
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Consequently, the degree of the second-order spatial coherence is 

9 {2) ( | Pi ~ P2 1 ) = 1 + somb 2 ( - - |pi - p 2 1 ) . (64) 

For a large value of 2R/d ~ A8, where A9 is the angular size of the radiation source viewed at the 
photodetectors, the point-spread somfe-function can be approximated as a (^-function of |pi — fel- 
We effectively have a "point-to-point" correlation between the transverse planes of Z\ = d and 
Z2 = d. In 1-D Eqs. (63 1 and (64 1 become 

G^( Xl - x 2 ) = It [1 + S mc 2 ( " A %^)] (65) 

A 



and 



g *Hxi-x*) = l + *r«' C Ae{Xl x - X2) ), (66) 



which has been experimentally demonstrated and reported in Fig. |18| 



We have thus derived the same second-order correlation and coherence functions as that of the 
quantum theory. The non-factorizable point-to-point correlation is expected at any intensity. The 
only requirement is a large number of point sub-sources with random relative phases participating 
to the measurement, such as trillions of independent atomic transitions. There is no surprise to 
derive the same result as that of the quantum theory from this simple model. Although the fields are 
not quantized and no quantum formula was used in the above calculation, this model has implied 
the same nonlocal two-photon interference mechanism as that of the quantum theory. Different 
from the phcnomenological theory of intensity fluctuations, this semiclassical model explores the 
physical cause of the phenomenon. 



5 Classical simulation 

There have been quite a few classical approaches to simulate type-one and type-two ghost 
imaging. Different from the natural non-factorizable type-one and type-two point-to-point imaging- 
forming correlations, classically simulated correlation functions are all factorizable. We briefly 
discuss two of these man-made factoriable classical correlations in the following. 

(I) Correlated laser beams. 

In 2002, Bennink et al. simulated ghost imaging by two correlated laser beams In this 
experiment, the authors intended to show that two correlated rotating laser beams can simulate 



the same physical effects as entangled states. Figure 21 is a schematic picture of the experiment 
of Bennink et al. Different from type-one and type-two ghost imaging, here the point-to-point 
correspondence between the object plane and the "image plane" is made artificially by two co- 
rotating laser beams "shot by shot" . The laser beams propagated in opposite directions and focused 
on the object and image planes, respectively. If laser beam-1 is blocked by the object mask there 
would be no joint-detection between D-y and D 2 for that "shot" , while if laser beam-1 is unblocked, 
a coincidence count will be recorded against that angular position of the co-rotating laser beams. 
A shadow of the object mask is then reconstructed in coincidences by the blocking— unblocking of 
laser beam-1. 

A man-made factorizable correlation of laser beam is not only different from the non-factorizable 
correlations in type-one and type-two ghost imaging, but also different from the standard statis- 
tical correlation of intensity fluctuations. Although the experiment of Bennink et al. obtained a 
ghost shadow, which may be useful for certain purposes, it is clear that the physics shown in their 
experiment is fundamentally different from that of ghost imaging. In fact, this experiment can be 
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Figure 21: A ghost shadow can be made in coincidences by "blocking-unblocking" of the 
correlated laser beams, or simply by "blocking-unblocking" two correlated gun shots. The 
man-made trivial "correlation" of either laser beams or gun shots are deterministic, i.e., the 
laser beams or the bullets know where to go in each shot, which are fundamentally different 
from the quantum mechanical nontrivial nondeterministic multi-particle correlation. 



considered as a good example to distinguish a man-made trivial deterministic classical intensity- 
intensity correlation from quantum entanglement and from a natural nonlocal nondeterministic 
multi-particle correlation. 

(II) Correlated speckles. 

Following a similar philosophy, Gatti et al. proposed a factorizable "speckle-speckle" classical 
correlation between two distant planes, p\ and P2, by imaging the speckles of the common light 
source onto the distant planes of p\ and P2 , [P3~] 

G (2) (/9i,p 2 ) oc 5(p - pi/m)6(p - ffa/m), (67) 

where po is the transverse coordinate in the plane of the light source)^] 

The schematic setup of the classical simulation of Gatti et al. is depicted in Fig. [22] [13 . Their 
experiment used cither entangled photon pairs of spontaneous parametric down-conversion (SPDC) 
or chaotic light for obtaining ghost shadows in coincidences. To distinguish from ghost imaging, 
Gatti et al. named their work "ghost imager" . The "ghost imager" comes from a man-made 
classical speckle-speckle correlation. The speckles observed on the object and image planes are the 
classical images of the speckles of the radiation source, reconstructed by the imaging lenses shown 
in the figure (the imaging lens may be part of a CCD camera used for the joint measurement). 
Each speckle on the source, such as the jth speckle near the top of the source, has two identical 
images on the object plane and on the image plane. Different from the non-factorizeable nonlocal 
image- forming correlation in type-one and type- two ghost imaging, mathematically, the speckle- 
speckle correlation is factorizeable into a product of two classical images of speckles. If two point 
photodetectors D\ and D2 are scanned on the object plane and the image plane, respectively, 
D\ and D2 will have more "coincidences" when they are in the position within the two identical 
speckles, such as the two jth speckles near the bottom of the object plane and the image plane. 
The blocking-unblocking of the speckles on the object plane by a mask will project a ghost shadow 
of the mask in the coincidences of £>i and D 2 . It is easy to see that the size of the identical speckles 

9 The original publications of Gatti et al. choose 2f-2f classical imaging systems with 1/2/ + 1/2/ = 1// to 
image the speckles of the source onto the object plane and the ghost image plane. The man-mde speckle-speckle 
image-forming correlation of Gatti et al. shown in Eq. |67| is factorizeable, which is fundamentally different from 
the natural non-factorizable imagc-formimg correlations in type-one and type-two ghost imaging. In fact, it is very 
easy to distinguish a classical simulation from ghost imaging by examining its experimental setup and operation. 
The man-made speckle-speckle correlation needs to have two sets of identical speckles observable (by the detectors 
or CCDs) on the object and the image planes. In thermal light ghost imaging, when using pseudo-thermal light 
source, the classical simulation requires a slow rotating ground grass in order to image the speckles of the source 
onto the object and image planes (typically, sub-Hertz to a few Hertz). However, to achieve a natural HBT non- 
factorizable correlation of chaotic light for type-two ghost imaging, we need to rotate the ground grass as fast as 
possible (typically, a few thousand Hertz, the higher the batter). 
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Figure 22: A ghost "imager" is made by blocking-unblocking the correlated speckles. The 
two identical sets of speckles on the object plane and the image plane, respectively, are the 
classical images of the speckles of the source plane. The lens, which may be part of a CCD 
camera used for the joint measurement, reconstructs classical images of the speckles of the 
source onto the object plane and the image plane, respectively. s and Si satisfy the Gaussian 
thin lens equation l/s a + 1/s, = 1//. 



determines the spatial resolution of the ghost shadow. This observation has been confirmed by 
quite a few experimental demonstrations. There is no surprise that Gatti et al. consider ghost 
imaging classical |27j . Their speckle-speckle correlation is a man-made classical correlation and 
their ghost imager is indeed classical. The classical simulation of Gatti et al. might be useful 
for certain applications, however, to claim the nature of ghost imaging in general as classical, 
perhaps, is too far |27j . The man-made factorizable speckle-speckle correlation of Gatti et al. is 
a classical simulation of the natural nonlocal point-to-point image-forming correlation of ghost 
imaging, despite the use of either entangled photon source or classical light. 



6 Local? Nonlocal? 



We have discussed the physics of both type-one and type-two ghost imaging. Although 
different radiation sources are used for different cases, these two types of experiments demonstrated 
a similar non-factorizable point-to-point image-forming correlation: 



Type-one: 



8{pi - p 2 ) 



dn s dPii S(k s + Hi) g\{K s , pi) g 2 (Ki, p 2 ) 



Type-two: 



i i 



1 2 

dndW \-^[g 1 (K,p 1 )g 2 (K ,p 2 ) + g 2 (K, p 2 )gx(tf ,pi)] \ ) 



(68) 



(69) 



Equations ( 68 ) and ( 69 1 indicate that the point-to-point correlation of ghost imaging, either type- 
one or type- two, is the results of two- photon interference. Unfortunately, neither of them is in the 
form of | J2j Ej\ 2 or \E\+E 2 1 2 , and neither is measured at a local space-time point. The interference 
shown in Eqs. (68) and (69 1 occurs at different space-time points through the measurements of 
two spatially separated independent photodetectors. 



In type-one ghost imaging, the <5-function in Eq. (68) means a typical EPR position-position 
correlation of an entangled photon pair. In EPR's language: when the pair is generated at the 
source the momentum and position of neither photon is determined, and neither photon-one nor 
photon-two "knows" where to go. However, if one of them is observed at a point at the object plane 
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the other one must be found at a unique point in the image plane. In type-two ghost imaging, 
although the position-position determination in Eq. ( 69 1 is only partial, it generates more surprises 
because of the chaotic nature of the radiation source. Photon-one and photon-two, emitted from a 
thermal source, are completely random and independent, i.e., both propagate freely to any direction 
and may arrive at any position in the object and image planes. Analogous to EPR's language: 
when the measured two photons were emitted from the thermal source, neither the momentum nor 
the position of any photon is determined. However, if one of them is observed at a point on the 
object plane the other one must have twice large probability to be found at a unique point in the 
image plane. Where does this partial correlation come from? If one insists on the view point of 
intensity fluctuation correlation, then, it is reasonable to ask why the intensities of the two light 
beams exhibit fluctuation correlations at pi = p 2 only? Recall that in the experiment of Sarcelli et 
al. the ghost image is measured in the near-field. Regardless of position, D\ and D 2 receive light 
from all (a large number) point sub-sources of the thermal source, and all sub-sources fluctuate 
randomly and independently. If A/1A/2 = for p\ ^ p 2 , what is the physics to cause A/1A/2 7^ 
at pi = p 2 ? 

The classical superposition is considered "local" . The Maxwell electromagnetic field theory 
requires the superposition of the electromagnetic fields, either | J2j Ej\ 2 or \Ei + E 2 \ 2 , takes place 
at a local space-time point (r,t). However, the superposition shown in Eqs. (68 1 and (69 1 happens 
at two different space-time points (ri,ti) and (r 2 ,t 2 ) and is measured by two independent pho- 
todetectors. Experimentally, it is not difficult to make the two photo-detection events space-like 
separated events. Following the definition given by EPR-Bell, we consider the superposition ap- 
pearing in Eqs. (68 1 and (69 1 nonlocal. Although the two-photon interference of thermal light can 
be written and calculated in terms of a semiclassical model, the nonlocal superposition appearing 
in Eq. (69) has no counterpart in the classical measurement theory of light, unless one forces a 
nonlocal classical theory by allowing the superposition to occur at a distance through the mea- 
surement of independent photodetectors, as we have done in Eq. (59 1. Perhaps, it would be more 
difficult to accept a nonlocal classical measurement theory of thermal light rather than to apply a 
quantum mechanical concept to "classical" thermal radiation. 



Conclusion: In summary, we may conclude that ghost imaging is the result of quantum in- 
terference. Either type-one or type-two, ghost imaging is characterized by a non-factorizable 
point-to-point image-forming correlation which is caused by constructive-destructive interferences 
involving the nonlocal superposition of two-photon amplitudes, a nonclassical entity corresponding 
to different yet indistinguishable alternative ways of producing a joint photo-detection event. The 
interference happens within a pair of photons and at two spatially separated coordinates. The 
multi-photon interference nature of ghost imaging determines its peculiar features: (1) it is non- 
local; (2) its imaging resolution differs from that of classical; and (3) the type-two ghost image 
is turbulence-free. Taking advantage of its quantum interference nature, a ghost imaging system 
may turn a local "bucket" sensor into a nonlocal imaging camera with classically unachievable 
imaging resolution. For instance, using the Sun as light source for type-two ghost imaging, we may 
achieve an imaging spatial resolution equivalent to that of a classical imaging system with a lens of 
92-meter diameter when taking pictures at 10 kilometers^] Furthermore, any phase disturbance 
in the optical path has no influence on the ghost image. To achieve these features the realization 
of multi-photon interference is necessary. 



Acknowledgment: The author thanks M. D Angelo, G. Scarcelli, J.M. Wen, T.B. Pittman, M.H. 
Rubin, and LA. Wu for helpful discussions. This work is partially supported by AFOSR and 
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10 The angular size of Sun is about 0.53°. To achieve a compatible image spatial resolution, a traditional camera 
must have a lens of 92-meter diameter when taking pictures at 10 kilometers. 
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Appendix: Presnel free-propagation 



We are interested in knowing how a known field E(ro,to) on the plane zq = propagates or 
diffracts into E(r,t) on another plane z — constant. We assume the field E(ro,to) is excited by 
an arbitrary source, either point-like or spatially extended. The observation plane of z = constant 
is located at an arbitrary distance from plane zq = 0, either far-field or near-field. Our goal is to 
find out a general solution E(r,t), or I(r,t), on the observation plane, based on our knowledge 
of E(ro,to) and the laws of the Maxwell electromagnetic wave theory. It is not easy to find such 
a general solution. However, the use of the Green's function or the field transfer function, which 
describes the propagation of each mode from the plane of zq = to the observation plane of 
z = constant, makes this goal formally achievable. 

Unless E(r , t ) is a non-analytic function in the space-time region of interest, there must exist 
a Fourier integral representation for _E(r ,to) 

E(r ,t ) = J dkE(k)w k (r ,t )e- wt0 , (A-l) 

where uik(i"o,io) is & solution of the Helmholtz wave equation under appropriate boundary con- 
ditions. The solution of the Maxwell wave equation u>k(ro,io) e~ luJt °, namely the Fourier mode, 
can be a s et of plane- waves or spherical- waves depending on the chosen boundary condition. In 
Eq. (A— 1|, E(k) = a(k)e 1(p ^ is the complex amplitude of the Fourier mode k. In principle we 



should be able to find an appropriate Green's function which propagates each mode under the 
Fourier integral point by point from the plane of zq = to the plane of observation, 

E(r,t) = J dkE(k)g(k : r-r ,t-to)w k (r Q ,to)e- tuJto 

= J dkg{k,r-r ,t-t )E{k,r ,t ), (A-2) 

where E(k,ro,to) = E(k) u>k(i"o, to) e~ luJt ° . The secondary wavelets that originated from each 
point on the plane of z = are then superposed coherently on each point on the observation plane 



with their after-propagation amplitudes and phases. It is convenient to write Eq. (A— 2) in the 
following form 

E(p,z,t) = j ' dujdKg(K,uj;p- po,z- z ,t - t ) E(k,uj; p , z ,t ), (^4-3) 

where we have used the transverse-longitudinal coordinates in space-time (p and z) and in mo- 
mentum (k, uj). 

Fig. |yl — l| is a simple example in which the field propagates freely from an aperture A of finite 



size on the plane oo to the observation plane a. Based on Fig. A—l we evaluate g(n, lo; p, z), namely 
the Green's function for free-space Fresnel propagation-diffraction. According to the Huygens- 
Fresnel principle the field at a given space-time point (p, z, t) is the result of a superposition of the 



spherical secondary wavelets that originated from each point on the <7q plane (see Fig. A—l ) 



E(p,z,t) = / diudn E(Z,w;Q,Q) I dp e -i(«t-*r) j (^_ 4 ) 



where we have set zq = and to — at plane do > aud defined r = y/ z 2 + \p — po\ 2 - In Eq. ( A— 4 1 , 
A(po) is the complex amplitude or relative distribution of the field on the plane of oo, which may 
be written as a simple aperture function in terms of the transverse coordinate po, as we have done 
in the earlier discussions. 

In the near-field Fresnel paraxial approximation, when \p — p \ 2 <C z 2 we take the first-order 
expansion of r in terms of z and p, 

r=^z 2 + \p- p \ 2 ~ z{\ + (A-5) 



a 



Figure A— I: Schematic of free-space Fresnel propagation. The complex amplitude A(p ) is 
composed of a real function A(po) and a phase e~ lK ' Po associated with each of the transverse 
wavevectors k in the plane of uq. Notice: only one mode of wavevector k(/?, lo) is shown in 
the figure. 



so that E(p,z,t) can be approximated as 



E(p,z,t)~ J dui dn E(k, w, 0, 0) J dpo — ^^e" c ~e" 

where e i 2^' /7 ~^ ' is named the Fresnel phase factor. 

Assuming that the complex amplitude A(po) is composed of a real function A(po) and a phase 
associated with the transverse wavevector and the transverse coordinate on the plane of 



do, as is reasonable for the setup of Fig. A— 1| we can then write E(p, z, t) in the form 
E(p,z,t) = f dujdn £;(K,w;0,0)e~ Mt — / dp A(p ) e^ e l 



The Green's function w; p, z) for free-space Fresnel propagation is thus 
g(jt,wj,z) = — [ dp A(p )e if{ -?°G(\p-p \,~). 



(A-6) 



In Eq. (A—6 \ we have defined a Gaussian function G(|a|,/3) = e 1 ^ 2 - 1 '"' , namely the Fresnel 
phase factor. It is straightforward to find that the Gaussian function G(|a|,/3) has the following 
properties: 

G*(\d\,P) = G(\a\,-p), 
GflSl.ft+A) = G{\d\,^)G(\a\^ 2 ), 
G{\d l +d 2 \,[3) = G{\d l \,f3)G{\d 2 \ 1 P)e^ s \ 



ddG(\d\,(3)e^ s = i^G{\%~). 



(A-7) 



Notice that the last equation in Eq. {A— 7) is the Fourier transform of the G(|a|,/3) function. As 
we shall sec in the following, these properties are very useful in simplifying the calculations of the 
Green's functions g(K,iv; p, z). 



Next, we consider inserting an imaginary plane a' between do and a. This is equivalent to 
having two consecutive Fresnel propagations with a diffraction-free a' plane of infinity. Thus, the 
calculation of these consecutive Fresnel propagations should yield the same Green's function as 
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that of the above direct Fresnel propagation shown in Eq. (^4—6): 



g(u,n;p,z) 

c 2 



d 1 d 2 



,!,,' j dp A(Po)G(\P - Pol, ^) G (\P- A ^) 



C 



dpoMpo)G(\p~ Po|,—) 



(AS) 



where C is a necessary normali zation constant for a valid Eq. (A— 8 I, and z = di + d 2 . The double 
integral of eipo and dp' in Eq. (A— 8) can be evaluated as 

dp' [ dp A(po) G(\p' - po|, G(\p - p'l 
Ja cdi cd 2 

dp A(po) G(p , G(p, ~) 
cdi cd 2 

dPG(p',-(^ + h)e-^ + ^' 
c d\ d 2 



i2nc d\d 2 
ui di+ d 2 



i2nc d\d 2 
lj di+ d 2 



dp A(Po) G(p , -j-) G(p, — ) 
. cdi cd 2 

di d 2 c d% + d 2 

dp A(Po)G(\p- po|, 73-— tt ) 
CT0 c («i + rf 2j 



where we have applied Eq. (A— 7 I, and the integral of dp' has been taken to infinity. Substituting 
this result into Eq. (A— 8) we obtain 



g(FZ,uj;p,z) 



C 2 



c 



i2nc e ^( dl+d2 ) 
lo di + d 2 



dp A(p )G(\p-po\, . 

a C («l + d 2) 



dpoA(p ) G(\p- po| 



cz 



Therefore, the normalization constant C must take the value of C — —iuj/2t:c. The normalized 
Green's function for free-space Fresnel propagation is thus 



g(n,uj;p,z) 



2nc z 



dpoA(po) G(\p- po|, — ). 

cz 



(A-9) 
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